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The Bethe-Salpeter equation restores exact elastic unitarity in the s— 
channel by summing up an infinite set of chiral loops. We use this equation 
to show how a chiral expansion can be undertaken in the two particle irre- 
ducible amplitude and the propagators accomplishing exact elastic unitarity 
at any step. Renormalizability of the amplitudes can be achieved by allow- 
ing for an infinite set of counter-terms as it is the case in ordinary Chiral 
Perturbation Theory. Crossing constraints can be imposed on the parame- 
ters to a given order. Within this framework, we calculate the leading and 
next-to-leading contributions to the elastic vrvr scattering amplitudes, for all 
isospin channels, and to the vector and scalar pion form factors in several 
renormalization schemes. A satisfactory description of amplitudes and form 
factors is obtained. In this latter case, Watson's theorem is automatically 
satisfied. From such studies we obtain a quite accurate determination of 
some of the ChPT 5C/(2)— low energy parameters (Ji — h = — 6.11q3 and 
Iq = 19.14 zb 0.19). We also compare the two loop piece of our amplitudes to 
recent two-loop calculations. 
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I. INTRODUCTION 



The dynamical origin of resonances in tttx scattering has been a recurrent subject in 
low energy particle physics Analyticity, unitarity, crossing and chiral symmetry have 
provided main insights into this subject. It is well known that so far no solution is avail- 
able exactly fulfilling all these requirements. In practical calculations some of the proper- 
ties mentioned above have to be given up. Standard Chiral Perturbation Theory (ChPT) 
furnishes exact crossing and restores unitarity order by order in the chiral expansion. In 
typical calculations the unitarity limit is reached at about a center of mass (CM) energy 
y/s ~ 4:\^f ~ 670 MeVQ still a low scale compared to the known resonances. Simply 
because of this reason standard ChPT is unable to describe the physically observed reso- 
nances, namely the p and the a. But even if the unitarity limit was much larger, and thus 
resonances would appear at significantly smaller scales than it, standard ChPT would be 
unable to generate them since its applicability requires the existence of a gap between the 
pion states and the hadronic states next in energy, which for tttt scattering are precisely 
the resonances. Resonances clearly indicate the presence of non perturbative physics, and 
thus a pole on the second Riemann sheet (signal of the resonance) cannot be obtained in 
perturbation theory to finite order. Hence, the energy regime for which ChPT holds has to 
satisfy s << m|. regardless of the unitarity limit. It so happens that both the resonances 
and the unitarity limit lead to similar scales, but there is so far no compelling reason to 
ascribe this coincidence to some underlying dynamical feature or symmetry]^. 

The desire to describe resonances using standard ChPT as a guide has led some authors 
to propose several approaches which favor some of the properties, that the exact scatter- 
ing amplitude should satisfy, respect to others. Thus, Pade Re-summation (PR) 0, Large 
A^/-Expansion (LNE) g]. Inverse Amplitude Method (lAM) [Q , Current Algebra Uni- 
tarization (CAU) 0, Dispersion Relations (DP) 0, Roy Equations ||^, Coupled Channel 
Lippmann-Schwinger Approach (CCLS) |@| and hybrid approaches have been suggested. 
Besides their advantages and success to describe the data in the low-lying resonance region, 
any of them has specific drawbacks. In all above approaches except by LNE and CCLS it 
is not clear which is the ChPT series of diagrams which has been summed up. This is not 
the case for the CCLS approach, but there a three momentum cut-off is introduced, hence 
breaking translational Lorentz invariance and therefore the scattering amplitude can be only 
evaluated in the the CM frame. On the other hand, though the LNE and CAU approaches 
preserve crossing symmetry, both of them violate unitarity. Likewise, those approaches 
which preserve exact unitarity violate crossing symmetry. 

A clear advantage of maintaining elastic unitarity lies in the unambiguous identification 
of the phase shifts. This is particularly useful to describe resonances since the modulus of 



Through this paper m is the pion mass, for which we take 139.57 MeV, and / the pion decay 
constant, for which we take 93.2 MeV. 

^ To see this, consider for instance the large Nc limit where / ~ \/iVc but ttir ~ A^^ and Tr ~ 1 /Nc, 
so in this limit / >> rriR, which suggests a possible scenario where resonances can appear well 
below the unitarity Umit. 
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the partial wave amplitude reaches at the resonant energy the maximum value allowed by 
unitarity. On the other hand, a traditional objection to any unitarization scheme is provided 
by the non-uniqueness of the procedure; this ambiguity is related to our lack of knowledge 
of an appropriate expansion parameter for resonant energy physics. This drawback does not 
invalidate the systematics and predictive power of unitarization methods, although it is true 
that improvement has a different meaning for different schemes. In addition, unitarization by 
itself is not sufficient to predict a resonance, some methods work while others do not. Against 
unitarization it is also argued that since crossing symmetry is violated, the connection to a 
Lagrangian framework is lost, and hence there is no predictive power for other processes in 
terms of a few phenomenological constants. 

The Bethe-Salpeter equation (BSE) provides a natural framework beyond perturbation 
theory to treat the relativistic two body problem from a Quantum Field Theory (QFT) point 



of view ||T0[ . This approach allows to treat both the study of the scattering and of the bound 
state properties of the system. In practical applications, however, approximations have to 
be introduced which, generally speaking, violate some known properties of the underlying 
QFT. This failure only reflects our in-capability of guessing the exact solution. The problem 
with the truncations is that almost always the micro-causality requirement is lost, and hence 
the local character of the theory. As a consequence, properties directly related to locality 
such as CPT, crossing and related ones are not exactly satisfied. This has led some authors 
to use the less stringent framework of relativistic quantum mechanics as a basis to formulate 



the few body problem |TT|. It would be, of course, of indubitable interest the formulation of 
a chiral expansion within such a framework. 

Despite of these problems, even at the lowest order approximation, or ladder approxi- 
mation, the BSE sums an infinite set of diagrams, allowing for a manifest implementation 
of elastic unitarity. This is certainly very relevant in the scattering region and more specif- 
ically if one aims to describe resonances, as we have discussed above. At higher orders, to 
implement unitarity one needs to include inelastic processes due to particle production ( 

7T7T —y KK, TTTT — >■ TTTTTTTr,- ■ ■). 

In general, the renormalization of the BSE is a difficult problem for QFT's in the con- 
tinuum. The complications arise because typical local features, i.e. properties depending 
on micro-causality, like crossing are broken by the approximate nature of the solution. This 
means that there is no a renormalized Lagrangian which exactly reproduces the amplitude, 
but only up to the approximate level of the solution. This is the price payed for manifest 
unitarity. Fortunately, from the point of view of the Effective Field Theory (EFT) idea, 
this problem can be tackled in a manageable and analytical way. Since this is a low energy 
expansion in terms of the appropriate relevant degrees of freedom, interactions are sup- 
pressed in powers of momentum and thus the BSE equation can be solved and renormalized 
explicitly by expanding the iterated two particle irreducible amplitude in a power series of 
momentum. This requires the introduction of a finite number of counter-terms for a given 
order in the expansion. The higher the order in this expansion, the larger the number of 
parameters, and thus the predictive power of the expansion diminishes. Thus an acceptable 
compromise between predictive power and degree of accuracy has to be reached. Within the 



2 



BSE such a program, although possible, has the unpleasant feature of arbitrariness in the 
renormalization scheme, although this is also the case in standard ChPT^ 

In this paper we study the scattering of pseudoscalar mesons by means of the BSE in the 
context of ChPT. Most useful information can be extracted from chiral symmetry, which 
dictates the energy dependence of the scattering amplitude as a power series expansion of 
1//^. The unknown coefficients of the expansion increase with the order. By using the BSE 
we can also predict the energy dependence of the scattering amplitude in terms of unknown 
coefficients and in agreement with the unitarity requirement. This can be done in a way to 
comply with the known behavior in ChPT. 

One important outcome of our calculation is the justification of several methods based on 
algebraic manipulations of the on-shell scattering amplitude, and which make no reference to 
the set of diagrams which are summed up. Moreover, some off-shell quantities are obtained. 
Although it is true that physical quantities only make sense when going to the mass shell 
there is no doubt that off-shell quantities do enter into few body calculations. 

Among others, the main results of the present investigation are the following ones: 

• A quantitative accurate description of both pion form factors and vrvr elastic scattering 
amplitudes is achieved in an energy region wider than the one in which ChPT works. 
We implement exact elastic unitarity whereas crossing symmetry is perturbatively 
restored. For the form factors, the approach presented here automatically satisfies 
Watson's theorem and it goes beyond the leading order Omnes representation which 
is traditionally used. 

• A meticulous error treatment is undertaken. Our investigations on the pion electric 
form factor together with our statistical and systematic error analysis provide a very 
accurate determination of some of the ChPT S'f/(2)— low energy parameters, 

- k = -G.ltoh h = 19.14 ± 0.19 



• The current approach reproduces ChPT to one loop. However and mainly due to the 
precise determination of the difference h — I2 quoted above, our one loop results have 
much smaller errors than most of the previously published ones. As a consequence, 
we generate some of the two-loop corrections more precisely than recent two loop 
calculations do. 

• The present framework allows for an improvement of any computed order of ChPT. 

The paper is organized as follows. In Sect. J| we discuss the BSE for the case of scat- 
tering of pseudoscalar mesons together with our particular conventions and definitions. In 
Sect. |Tl|we review and improve the lowest order solution for the off-shell amplitudes, found 



in a previous work |12]. We also show that, already at this level, a satisfactory description 



^There, higher order terms than the computed ones are exactly set to zero. This choice is as arbi- 
trary as any other one where these higher order terms are set different from zero. The advantages 
of taking them different from zero as dictated by unitarity will become clear along the paper. 
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of data can be achieved when a reasonable set of ChPT low energy parameters is used. This 
is done in Subsect. |III F| . Technical details are postponed to Appendix ^ The computation 
of the next-to-leading order corrections, within this off-shell scheme, turns out to be rather 
cumbersome. The difficulties can be circumvented by introducing what we call the "on- 
shell" scheme, as we do in Sect. 0. After having discussed some renormalization issues, the 
decoupling of the off- shell amplitudes is achieved in Subsect. [IV C| by means of a generalized 
partial wave expansion. In this context, the renormalization conditions can also be stated in 
a more transparent way. The findings of Subsect. |IV C| allow us to write an exact T— matrix 
if an exact two particle irreducible amplitude was known. In Subsect. |1V D| , we show that the 
found solution satisfies off-shell unitarity and in the next two subsections we present several 
systematic expansions of the two particle irreducible amplitude. In Subsect. |IV G| , we show 
how our amplitudes can be fruitfully employed for the calculation of pseudoscalar meson 
form factors, in harmony with Watson's theorem. Next-to-leading order numerical results 
for both TTTT scattering and form factors (vector and scalar) are presented in Subsect. IV H . 
This study allows us to determine, very precisely, the parameters hyh (specially their dif- 
ference) and Iq from experimental data. Besides, we also compare with recent two-loop 
calculations. In Subsect. |IV1| we compare some results obtained at leading and next-to- 
leading accuracy, within the BSE on-shell scheme. Such a study tests the convergence of 
the approach presented in this work. In Sect. we try to understand qualitatively the origin 
of some renormalization constants which naturally appear within the BSE approach and in 
particular a remarkable formula for the width of the p meson, very similar to the celebrated 
KSFR, is deduced. In Sect. ^ we briefly analyze alternative unitarization methods on the 
light of the BSE approach. Conclusions are presented in Sect. |VI1] . Finally, in Appendix ^ 
the leading and next-to-leading elastic tttt scattering amplitudes in ChPT are compiled and 
we also give analytical expressions for the u— and t— unitarity chiral corrections projected 
over both isospin and angular momentum. 



II. THE BETHE-SALPETER EQUATION. 

Let us consider the scattering of two identical mesons of mass, m. The BSE can, with 
the kinematics described in Fig. |l], be written as 

Tp{p, k) = Vp{p, k) + ij -0j-Jp{q, A;) A(g+) A(g_) V>(p, q) (1) 

where q± = (-P/2 ± q) and Tp{p,k) and Vp{p,k) are the total scattering amplitude^ and 
the two particle irreducible amplitude (potential) respectively. Besides, A is the exact pseu- 
doscalar meson propagator. Note that the previous equation requires knowledge about the 



^The normalization of the amplitude T is determined by its relation with the differential cross 
section in the CM system of the two identical mesons and it is given by da/dCl = \Tp{p, A;)p/647r^s, 
where s = P^. The phase of the amplitude T is such that the optical theorem reads lmTp{p,p) = 
— (Ttot('S^ — 4s m^)^/^, with fitot the total cross section. The contribution to the amputated Feynman 
diagram is (— iTp(p, k)) in the Fig. ||. 
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off-shell potential and the off-shell amplitude^. Clearly, for the exact potential V and propa- 
gator A the BSE provides an exact solution of the scattering amplitude T Obviously an 
exact solution for T is not accessible, since V and A are not exactly known. An interesting 
property, direct consequence of the two particle irreducible character of the potential, is that 
in the elastic scattering region s > 4m^, is a real function, and that it also has a discon- 
tinuity for s < 0, i.e. t > Am?. We will also see below that, because of the inherent freedom 
of the renormalization program of an EFT, the definition of the potential is ambiguous, and 
some reference scale ought to be introduced in general. The exact amplitude, of course, will 
be scale independent. 




FIG. 1. Diagrammatic representation of the BSE equation. It is also sketched the used kinematics. 

Probably, the most appealing feature of the BSE is that it accomplishes the two particle 
unitarity requirement^ 

Tp(p, k) - Tp{k,pr = -i(2vr)2 J ^Jp{q. k) 5+ {ql - m') 6+ {ql - m') Tpiq^pf (2) 

as can be deduced from the real character of the potential above threshold s > 4m^, with 
6~^{p'^ — im?) = Q{p^)6{p'^ — m?). It is important to stress that although the above unitarity 
condition is considered most frequently for the on-shell amplitude, it is also fulfilled even for 
off-shell amplitudes. We will see below that these conditions are verified in practice by our 
amplitudes. 

Isospin invariance leads to the following decomposition of the two identical isovector 
meson scattering amplitude for the process (-P/2 + p,a) + (-P/2 — p,b) — >• (-P/2 + k,c) + 
{P/2-k,d), 

Tp{p, k)ab-cd = Ap{p, k)5ab5cd + Bp{p, k)5achd + Cp{p, k)5adhc (3) 

with a, 6, c and d Cartesian isospin indices. With our conventions, the Mandelstam variables 
are defined as s = P^, t = {p — k)"^ and finally u = {p + k^ and the isospin projection 
operators, P^, are given by [|l^] 



^We will see later that this off-shellness can be disregarded within the framework of EFT's and 
using an appropriate regularization scheme. 

^Cutkosky's rules lead to the substitution A(p) — >= (— 27ri)J"*"(p^ — rm?) 
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pO — A 

^ab;cd — r^'Jab<Jcd 

Plb;cd = i^acSbd - SadSbc) 

Pab;cd = (^^aJbd + ^adhc " ^^^ab^cd^ (4) 

The isospin decomposition of the amphtude is then Tp(p, k)ab;cd = J2i Plb-cd'^piP^ where 

TUp, k) = l {3Apip, k) + Bp{p, k) + Cp{p, k)) 
T},{p,k) = ^iBpip,k)-Cpip,k)) 

TUp,k) = ^{Bp{p,k) + Cp{p,k)) (5) 

Based on the identical particle character of the scattered particles we have the following 
symmetry properties for the isospin amplitudes 

n{p, k) = {-ytUp, -k) = i-ym-p, k) (6) 

The above relations imply: 

Ap{p, k) = Ap{p, -k) = Ap{-p, k) 

Bp{p, k) = Cp{p, -k) = Cp{-p, k) (7) 
On the other hand crossing symmetry requires 

rr ( i.\ ^ f p + k + P p + k-P\ 

Tp{p, k)ab;cd = T^p-k) I ^ , ^ I (8) 

\ / ac;bd 

which implies 

fp + k + P p + k- P\ 



Bp{p, k)=A 



(p-k) 



2 ' 2 
'p + k + P p + k- P" 



Cp{p, k) = ^ , '—^ j (9) 

Finally crossing, together with rotational invariance also implies: 

Tp{p, k)ab-cd = T_p{-k, -p)cd;ab = T_p{-k, -p)ab;cd = Tp{k, p)ab;cd (lO) 

which forces to all three functions Ap,Bp and Cp to be symmetric under the exchange 
of k and p. The relations of Eqs. (|^) and can be combined to obtain the standard 
parameterization on the mass shell: 

Ap{p, k) = A{s, t, u) = A{s, u, t) 
Bp{p,k) =A{t,s,u) 

Cp{p,k)=A{u,t,s) (11) 
6 



On the other hand the relations in Eq. can be inverted, and thus one gets 

Ap{p,k) = '^{TUp,k)-TUpJ 

Bp{p,k)=TUp,k)+T}.{p,k) 
Cp{p,k)=T'p{p,k)-TUp,k) 



(12) 



Hence, the crossing conditions stated in Eq. ([lID , impose a set of non trivial relations between 
the amplitudes Tp{p, k), Tp{p, k) and Tp(p, k). 

For comparison with the experimental CM phase shifts, we define the on-shell 

amplitude for each isospin channel as 



T\s,t)=Tl{p,k) 



P 



m 



s/4; 



P .p = P . k = 



(13) 



and then the projection over each partial wave J in the CM frame, Tjj{s), is given by 



TiAs) 



1 r+i 



2 J- 



Pj {cos 6) Tp{p, k) d{cosl 



Aa (s, m^, m^) 



.2i<5/j(s) 



(14) 



where 6 is the angle between p and k in the CM frame, Pj the Legendre polynomials and 
X{x, y, z) = + + — 2xy — 2xz — 2yz. Notice that in our normalization the unitarity 
limit implies \Tjj{s)\ < IGns/X^ 



III. OFF-SHELL BSE SCHEME 

A. Chiral expansion of the potential and propagator 

We propose an expansion along the lines of ChPT both for the exact potential {V) and 
the exact propagator (A), 

A(p) = A(°)(p) + A(2)(p) + ... 
Vp{p, k) = Vp(p, k) +(2) Vp(p, k) + ... (15) 

Thus, at lowest order in this expansion, V should be replaced by the 0{p^) chiral amplitude 
((2)T) and A by the free meson propagator, A'^°)(r) = (r^ — m^ + ie) ^. Even at lowest order, 
by solving Eq. (P we sum up an infinite set of diagrams. This approach, at lowest order 
and in the chiral limit reproduces the bubble re-summation undertaken in Ref. |T^. This 
expansion is related to the approach recently pursued for low energy iViV— scattering where 
higher order t— and m— channel contributions to the potential are suppressed in the heavy 
nucleon mass limit [ 1^ . 



To illustrate the procedure, let us consider elastic vrvr scattering in the 1 = 0,1,2 channels. 
One word of caution should be said, before proceeding to practical calculations. As one sees 
the solution of the BSE requires knowledge about off-shell quantities, like the potential, 
Vp{p^ k) . This is generally an ambiguous quantity which depends on the particular choice 
of the pion field. On- shell quantities should however be independent of such a choice. The 
calculations carried out below correspond to a particular ansatz. A more detailed discussion 
on other possible choices is postponed to Appendix [A6| . The main result is that a judicious 
choice of the renormalization scheme, leads to on-shell quantities free from this arbitrariness. 
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1. 1=0 TTTT scattering. 



At lowest order, the off-shell potential V in this channel is given by 

To solve Eq. (|I]) with the above potential we proceed by iteration. The second Born approx- 
imation suggests a solution of the form 

T°(p, k) = A{s) + B{s){p' + + C{s)p^k'' (17) 

where A[s),B[s) and C{s) are functions to be determined. Note that, as a simple one loop 
calculation shows, there appears a new off-shell dependence (p^fc^) not present in the 0{p^) 
potential ^^-"T. That is similar to what happens in standard ChPT 

The above ansatz reduces the BSE to a linear algebraic system of four equations with 
three unknowns (Eq. ( [Al| ) in Appendix ^). The system turns out to be compatible and the 
solution of it is given in Eqs. ( [A5| ) and (|A6|) of the Appendix 0. 



At the lowest order in the chiral expansion examined here, the isoscalar amplitude on 
the mass shell and in the CM frame (P = 0,p° = A;*^ = 0, P° = ^/s) is purely s— wave, and 
its inverse can be obtained from Eqs. ( [A5| ) and ( |A6| ) . It reads 

00 KS) olsj + 2h{Aw?) + {w? - 2s)p + {s - Am^) hiAm^) ^ ^ 

where Io{s) is a logarithmically divergent integral, which explicit expression is given in 
Eq. ( |A2| ) . Similarly /2(4m^) and /4(4m^) are divergent quantities, which are defined in 
terms of the quadratic and quartic divergent integrals his) and h{s) also introduced in 
Eq. (|A^). Thus the above expression for the inverse amplitude requires renormalization, we 
will back to this point in Subsect. ^11 D| . 



2. 1=1 TTTT scattering 
At lowest order, the off-shell potential V in this channel is approximated by 

vhip^k) nip^k) = ^-^ (19) 

As before, to solve Eq. (||) with the above potential we propose a solution of the form 

r>(p, k) = M{s)p ■ k + N{s){p ■P){k- P) (20) 

where M and N are functions to be determined. Note that, as expected from our previous 
discussion for the isoscalar case, there appears a new off-shell dependence {{p ■ P){k ■ P)) 
not present in the 0{p'^) potential. Again, this ansatz reduces the BSE to a linear algebraic 
system of equations which provides the full off-shell scattering amplitude, which is given in 
Eq. (^). 
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At the lowest order presented here, we have only j9— wave contribution and the resulting 
inverse CM amplitude on the mass shell, after angular momentum projection, reads 



Tu\s) = = -Io{s) + (21) 
Similarly to the isoscalar case discussed previously, the above equation presents divergences 



which need to be consistently renormalized, this issue will be addressed in Subsect. Ill D. 



3. 1=2 vrvr scattering. 
At lowest order, the off-shell potential V in this channel is given by 

k) -(2) TUp, k) = " ^f, + (22) 

This resembles very much the potential for the isoscalar case, and we search for a solution 
of the BSE of the form 

T|(p, k) = A{s) + B{s)ip' + + C{s)p'e (23) 

Similarly to the case / = 0, the functions A{s), B{s) and C{s) can be readily determined 
and are given in the Appendix 0. 

At the lowest order in the chiral expansion examined here, the 1 = 2 amplitude on the 
mass shell and in the CM frame is purely s— wave, and from Eqs. ( |A9|) and ( [AlOp ) we find 
its inverse reads 

T-Us) = -I (s) + 2iP + U4m^)f 

20 y^) 2/4(4m2) + (s-2m2)/2 + (s-4m2)/2(4m2) ^ ^ 

Once again, the above equation has to be renormalized. 



B. On-shell and off-shell unitarity 

As we have already anticipated, the solutions of the BSE must satisfy on-shell and off- 
shell unitarity. This is an important check for our amplitudes. This implies in turn conditions 
on the discontinuity (Disc [/(s)] = f{s + ie) — f{s — ie), s > Aw?) of the functions v4(s), 
B[s) and C{s) for the / = and 1 = 2 cases and N{s) and M{s) for the isovector one. 
Going through the unitarity conditions implicit in Eq. (0), we get a set of constraints which 
are compiled in Eqs. ([All|) and (|A12|) of the Appendix 0. 



After a little of algebra, one can readily check that the discontinuity conditions of Eqs. ( |A11 



( |A12|) are satisfied by the off-shell amphtudes found in Subsect. [Ill A| . That guarantees that 



the solutions of the BSE found in that subsection satisfy both off-shell and on-shell unitarity. 

For the on-shell case, elastic unitarity can be checked in a much simpler manner than that 
presented up to now. The on-shell amplitudes can be expressed in the following suggestive 
form which, as we will see, can be understood in terms of dispersion relations. 



9 



Tfiis) = -lois) - Cjj + —— (25) 

yij[s) 

where Cjj is a constant and the potentials are trivially read off from the on shell amplitudes 
in Eqs (|TBp,(|2T]) and (P^. These potentials contain an infinite power series of 1//^. In the 
on-shell limit, the unitarity condition for the partial waves is more easily expressed in terms 
of the inverse amplitude, for which the optical theorem for s > 4m^ reads 



ImT-(. + ie) = = -^./l-l^ = -Im/o(. + ie) (26) 



thus, the on-shell amplitudes found in Subsect. [Ill A| for the several isospin- angular momen- 
tum channels trivially meet this requirement. 



C. Crossing properties 

On the mass-shell and at the lowest order in the chiral expansion presented in Sub- 
sect. ^TM the Eq. (|12D leads to 



2 

Ap{p, k) = -(Too(s) - T2o{s] 

u — t 

Bp{p, k) = T,ois) - 3 -^Tnis) 

s — 4m^ 

V — f 

Cpip, k) = T2ois) + 3j-^Tu{s) (27) 



Obviously, the crossing conditions stated in Eq. (|TT]) are not satisfied by the functions defined 
in Eq. ( P?] ) , although they are fulfilled at lowest order in 1//^ . This is a common problem 
in all unitarization schemes. In the BSE the origin lies in the fact that the kernel of the 
equation, / d'^q ■ ■ ■ A(g+) A(g_) ■ ■ ■ breaks explicitly crossing, and given a potential V it only 
sums up all s— channel loop contributions generated by it. Thus this symmetry is only 
recovered when an exact potential V, containing t— and m— channel loop contributions, is 
used. 

Actually, in our treatment of the channels / = 0, 1, 2 we require 3+2+3 = 8 undetermined 
constants^ , whereas crossing imposes to order 1//^ only 4, namely ^i, 2,3,4- This is not as 
severe as one might think since we are summing up an infinite series in 1//^. Nevertheless, 
we will show below how this information on crossing can be implemented at the level of 
partial waves. 



we will discuss in the next subsection, a consistent renormalization program allows one to 
take the divergent integrals Io{4:m?), I2{'im^), l4:{4:m?) independent of each other in each isospin 
channel. 
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D. Renormalization of the amplitudes 



To renormalize the on-shell amplitudes given in Eqs. (P^) , (pi]) and (P^, or the cor- 
responding ones for the case of off-shell scattering, we note that in the spirit of an EFT 
all possible counter-terms should be considered. This can be achieved in our case in a 
perturbative manner, making use of the formal expansion of the bare amplitude 

T = V + VGqV + VGqVGqV + ■ ■ ■ , (28) 

where Gq is the two particle propagator. Thus, a counter-term series should be added to the 
bare amplitude such that the sum of both becomes finite. At each order in the perturbative 
expansion, the divergent part of the counter-term series is completely determined. However, 
the finite piece remains arbitrary. Our renormalization scheme is such that the renormalized 



amplitude can be cast, again, as in Eqs. (|Tq),(pTl) and (P4|) . This amounts in practice, 
to interpret the previously divergent quantities /2n(4m^) as renormalized free parameters. 
After having renormalized, we add a superscript R to differentiate between the previously 
divergent, J2n(4m^), and now finite quantities, J|^(4m^). These parameters and therefore the 
renormalized amplitude can be expressed in terms of physical (measurable) magnitudes. In 
principle, these quantities should be understood in terms of the underlying QCD dynamics, 
but in practice it seems more convenient so far to fit /^(4m^) to the available data. The 
threshold properties of the amplitude (scattering length, effective range, etc..) can then be 
determined from them. Besides the pion properties m and /, at this order in the expansion 
we have 8 parameters. The appearance of 8 new parameters is not surprising because the 
highest divergence we find is quartic {Ii{s)) for the channels J = and 1 = 2 and quadratic 
(12(5)) for the isovector channel and therefore to make the amplitudes convergent we need 
to perform 3+3-1-2 subtractions respectively. This situation is similar to what happens in 
standard ChPT where one needs to include low-energy parameters (Z's). In fact, if t— and 
u— channel unitarity corrections are neglected, a comparison of our (now) finite amplitude, 
Eqs. (p^,(|2TD and (p^, to the O(p^) vrvr amplitude in terms of some of these Ps becomes 
possible. Such a comparison will be discussed in the next subsection. 

This renormalization scheme leads to a renormalized amplitude which does not derive 
from a renormalized Lagrangian. This is again closely linked to the violation of exact crossing 
symmetry, and it is detailed discussed in the Sect. |A 5| of the Appendix |^. 



E. Crossing Symmetry Restoration and Comparison with one loop ChPT 

Our amplitudes contain undetermined parameters which, as stated previously outnumber 
those allowed by crossing symmetry at 0{p'^). Nevertheless, by imposing suitable constraints 
on our parameters we can fulfill crossing symmetry approximately. We do this at the level 



of partial wave amplitudes. For completeness, we reproduce here a discussion from Ref. [12 
where this issue was first addressed. At the lowest order in the chiral expansion proposed 
in this section, we approximate, in the scattering region s > 4m^, the (9(1//"^) t— and u— 
channel unitarity corrections (function hjj in Eq. (p^) of the Appendix |BD by a Taylor 
expansion around threshold to order [s — 4m^)^. At next order in our expansion (when 
the full 0(p'') -corrections are included both in the potential and in the pion propagator) 
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we will recover the full t— and u— channel unitarity logs at 0{l/f^), and at the next 
order (0(1//^)), we will be approximating these logs by a Taylor expansion to order (s — 
4m^)^. Thus, the analytical structure of the amplitude derived from the left hand cut is 
only recovered perturbatively. This is in common to other approaches (PR, DP, lAM ■ ■ •) 
fulfilling exact unitarity in the s-channel, as discussed in [^, 0, ||T^. Thus, this approach 
violates crossing symmetry. At order (9(1//'^) our isoscalar s— , isovector p— and isotensor 
s— wave amplitudes are polynomials of degree two in the variable {s — 4m^), with a total 
of eight (3+2+3) arbitrary coefficients {Iq^2,a (4m^), /o,2 ~ (4m^), ^0,2,4^ (4m^) ), and there 
are no logarithmic corrections to account for t— and m— channel unitarity corrections. Far 
from the left hand cut, these latter corrections can be expanded in a Taylor series to order 
(s— 4m^)^, but in that case the one loop SU{2) ChPT amplitudes can be cast as second order 
polynomials in the variable (s — Aim?), with a total of four (/i,2,3,4) arbitrary coefficients |]16 



To restore, in this approximation, crossing symmetry in our amplitudes requires the existence 
of four constraints between our eight undetermined parameters. These relations can be found 
in Eq. ( |A15|) in the Appendix ^ 

Once these constraints are implemented in our model, there exists a linear relation be- 
tween our remaining four undetermined parameters (/(f/~°(4m^), J(f2^~^(4m^)) and the most 
commonly used /i, 2,3,4 parameters. Thus, all eight parameters (J(^2fr°'^(4'^^), -^(f/~^(4m^)) 
can be expressed in terms of Ti, 2,3,4 (see Eq. ( |A16| ) of Appendix 0). 



F. Numerical results for / = 0, 1, 2 



After the above discussion, it is clear that at this order we have four independent pa- 
rameters /(f'^~^ and l2'^~^ which can be determined either from a combined 
X^— fit to the isoscalar and isotensor s— and isovector p— wave elastic tttt phase shifts or 
through, Eq. ([A16| ), from the Gasser-Leutwyler or other estimates of the ^i, 2,3,4 low energy 



parameters. In a previous work, ||T2[, we have already discussed the first procedure (x^— fit) 
and thus we will follow here the second one. Therefore, we will try to address the follow- 
ing question: Does the lowest order of the off-shell BSE approach together with reasonable 
values for the ^i, 2,3,4 parameters describe the observed vrvr phase-shifts in the intermediate 
energy region? To answer the question, we will consider two sets of parameters: 

set A: h = -0.62 ± 0.94, k = 6.28 ± 0.48, [3 = 2.9 ± 2.4, k = 4:A ± 0.3 

set B: h = -1.7 ± 1.0 , k = 6.1 ± 0.5 , [3 = 2.9 ± 2.4, h = 4:A ± 0.3 (29) 

In both sets I3 and Z4 have been determined from the 5*^7(3) mass formulae and the scalar 



radius as suggested in [|16[ and in [|T8|, respectively. On the other hand the values of li^2 come 



from the analysis of Ref. [IS] of the data on K14— decays (set A) and from the combined study 
of K14— decays and vrvr with some unitarization procedure (set B) performed in Ref. pO| . 

Results are presented in Fig. ^ and in Table |. In the figure we show the prediction (solid 
lines) of the off-shell BSE approach, at lowest order, for the s— and wave vrvr scattering 
phase-shifts for all isospin channels and for both sets A (left panels) and B (right panels) of 
the Us parameters. We assume Gauss distributed errors for the Z's parameters and propagate 
those to the scattering phase shifts, effective range parameters, etc. . . by means of a Monte 
Carlo simulation. Central values for the phase-shifts have been computed using the central 
values of the Vs parameters. Dashed lines in the plots are the 68% confidence limits. 



12 



As we see for both sets of constants, the simple approach presented here describes the 
isovector and isotensor channels up to energies above 1 GeV, whereas the isoscalar channel 
is well reproduced up to 0.8-0.9 GeV. In the latter case, and for these high energies, one 
should also include the mixing with the KK channel as pointed out recently in Refs. 0. 
Regarding the deduced threshold parameters we find agreement with the measured values 
when both theoretical and experimental uncertainties are taken into account. Furthermore, 
both sets of parameters predict the existence of the p resonanc^ in good agreement with 
the experimental data, 

set A: mp = TTOl^g [MeV], Vp = 180+1^ [MeV] 
set B: mp = 715tlo [MeV], Tp = 130^1^ [MeV] 

(31) 

In this way the "existence" of the p resonance can be regarded as a prediction of the BSE 
with ChPT and the parameters obtained from some low energy data. Favoring one of the 
considered sets of parameters, implies an assumption on the size of the 0{p^) contributions 
not included at this level of approximation. Finally, it is worth mentioning that the pre- 
dicted parameters J^(4m^), agree reasonably well with those fitted to tttt scattering data in 
Ref. [0. 





J = 
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J = 


/ = 1 


J = 0,/ 


= 2 




set A 


set B 


set A 


set B 


set A 


set B 


-102lo«(4m2) 


3.1(4) 


2.6(5) 


9.8(14) 


10.9(14) 


6.6(6) 


6.1(6) 


-103l^(4m2) 


1.7(11) 


3.0(12) 


-77(7) 


-83(8) 


7.5(12) 


8.4(12) 


-lO^/f (W) 


3.0(12) 


2.7(12) 






2.1(5) 


2.1(5) 




225(7) 


218(7) 


39.5(11) 


39.5(12) 


-41.0(12) 


-42.1(13) 




308(20) 


286(20) 


7.3(14) 


8.6(16) 


-72.5(20) 


-74.4(23) 



TABLE I. OfT-shcll BSE approach parameters (/^(4to^)) obtained from both sets of Ts parameters given 
in Eq. (p9[). I^{A'm?) are given in units of (2171)". We also give the threshold parameters ajj and bjj obtained 
from an expansion of the scattering amplitude KeTjj = — 167rm(s/4 — m^)'^[a/j + bij{s/A — m?) + • • •] 
close to threshold. Errors have been propagated by means of a Monte Carlo simulation, they are given in 
brackets and affect to the last digit of the quoted quantities. 



^We determine the position of the resonance, by assuming zero background and thus demanding 
the phase-shift, 5ii, to be 7r/2. Furthermore, we obtain the width of the resonance from 



1 



{m?p — s) tan 5ii(s) 



nir. 



d5ii{s) 



ds 



IQ-Km^p 



dRerf/(s) 



\'i{rn'jj,m? ,m?) 



ds 



(30) 
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FIG. 2. Several tttt phase shifts as a function of the total CM energy for both sets of /'s quoted 
in Eq. (p9|). Left (right) panels have been obtained with the set A (B) of parameters. Solid lines are 
the predictions of the off-shell BSE approach, at lowest order, for the different U— channels. Dashed lines 
are the 68% confidence limits. Top panels (/ = 0, J = 0): circles stand for the experimental analysis of 
Refs. - ||2^. Middle panels (/ = 1, J = 1): circles stand for the experimental analysis of Refs. |2l| and 
p3| . Bottom panels {I — 2, J — 0): circles stand for the experimental analysis of Ref. In all plots 

the triangles are the Frogatt and Petersen phase-shifts (Ref. |2^]) with no errors due to the lack of error 
estimates in the original analysis. 
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IV. ON-SHELL BSE SCHEME 



As already recognized in previous work [|T^, it is very hard in practice to pursue the 



calculation of higher orders, within the off-shell scheme presented in the previous section. 
The difficulty lies on the non-locality of the four point Green-functions involved in the BSE in 
momentum space, and the subsequent and unavoidable renormalization. More specifically, 
to keep the full off-shell dependence of V and T in the BSE is the origin of most of the 
difficulties. 



However, if we look at the amplitude in each isospin channel separately (Eqs. (p^,(|2TD 
and (p^ ) it seems that the off-shellness can be ignored by simply renormalizing the pa- 
rameters of the lowest order amplitude (/ and m) . The authors of Refs. ^ take advantage 
of this observation. Though, it is important to stress here, that this renormalization is 
different for each JJ— channel and therefore this procedure does not provide a satisfactory 
renormalization scheme. 

In this section, we come up with a consistent and computationally feasible renormaliza- 
tion program where all off-shellness effects can be incorporated into a legitimate redefinition 
of the two particle irreducible amplitude, V, at all orders of the chiral expansion. 

This new scheme allows us not only to describe the data but also to make predictions 
for some of the ChPT-low energy constants. 



A. Renormalization of power divergences 

Let us consider the iteration in the BSE of a renormalized potential, instead of the "bare" 
potential . Obviously, to define such a quantity requires a renormalization prescription]^. 
This new viewpoint can only be taken, as we will see, because the divergence structure of 
the non-linear sigma model is such that, all non- logarithmic divergences can be absorbed 
by a suitable renormalization of the parameters of the Lagrangian. For completeness, we 
repeat here the argument given in Ref. in a way that it can be easily applied to our case 



of interest. It is easier to consider first the chiral limit m = 0. The Lagrangian at lowest 
order is given by 



2 



C = Ltr(^d''U^d^,u) (32) 

with U a dimensionless unitarity matrix, involving the Goldstone fields, which transforms 
linearly under the chiral group. Thus, the necessary counter-terms to cancel L loops are 
suppressed by the power /^Ci--^). Let A be a chiral invariant regulator with dimension of 
energy and D the dimension of a counter-term which appears at L loops. Since the U field 
is dimensionless, D simply counts the number of derivatives. The operator appearing in the 
counter-term Lagrangian should have dimension 4, thus we have 

D = 2 + 2L-r (33) 



^In the following we will assume a mass independent regularization scheme, such as e.g. dimen- 
sional regularization. The reason for doing this is that it preserves chiral symmetry. 
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with r the number of powers of the regulator A which accompany this counter-term in the 
Lagrangian (degree of divergence). If we denote the counter-terms in the schematic way 

oo 2L 

c.. = E f^'~'' E * E A'{d'^''-n^ (34) 

L=0 i r=0 

where {d^)i denotes a set of chiral invariant hnear independent operators made out of the 
matrix U and comprising D derivatives and cf^ are suitable dimensionless coefficients. We 
see that at L loops the only new structures are those corresponding to A° (r = 0) which 
in actual calculations corresponds to logarithmic divergences. Higher powers in A generate 
structures which were already present at L — 1 loops. For instance, let us consider the 
one-loop correction, by construction it is down respect to the leading order by two powers 
of /, i.e., the corresponding counter-terms in the Lagrangian will be made out of chiral 
invariant operators constructed by means of derivatives of the dimensionless U matrix, 
defined above, and any of them will be multiplied by an overall factor /°. In addition, 
the total energy dimension of the counter-terms has to be four. Hence, only terms with two 
and four derivatives[^ might appear. To get the proper dimensions, these terms will have 
to be multiplied by A^ and A° In A respectively. Thus, the quadratic divergence, first term, 
is a renormalization of the leading kinetic energy piece (Eq. (|32D). Besides, to renormalize 
the logarithmic divergence, one has to add new counter-terms, not present in the original 
Lagrangian (Eq. (p2D). For higher loops, the argument can be easily generalized, see for 
more details Ref. |]30| . 

This result is very important because it means that all but the logarithmic divergences 
can effectively be ignored. An efficient scheme which accomplishes this property is dimen- 
sional regularization, since 




This argument can be extended away from the chiral limit, although in this case new terms 
appear which vanish in the limit 0. The conclusion is again the same. If every possible 

counter-term compatible with chiral symmetry is written down, all but the logarithmic 
divergent pieces can be ignored. 

The above discussion means that when renormalizing the BSE, we may set to zero (at a 
given renormalization point) all power divergent integrals which appeared in Sect. This 
is because they only amount to a renormalization of the undetermined parameters of the 
higher order terms of the Lagrangian. The power of this result is a direct consequence of 
the symmetry, the derivative coupling of the pion interactions and the fact that in ChPT 
there is an infinite tower of operators. In a sense this is only true for the exact theory 
with an infinite set of counter-terms, it is to say when we iterate by means of the BSE a 
"renormalized" two particle irreducible amplitude, l^, with an infinite number of terms. 



Note that an odd-number of derivatives are forbidden by parity conservation and the only 
chirally invariant operator with zero derivatives is tr(C/{7^) = 2, an irrelevant constant. 
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To keep all power divergences and not only the logarithmic ones and simultaneously 
iterate the most general two particle irreducible amplitude leads to redundant combina- 
tions of undetermined parameters. For instance, to include in the potential the tree 0{p^) 
Lagrangian (Ps) and keep the power (non logarithmic) divergences produced by the loops 
made out of the 0{p'^) pieces of the Lagrangian produces redundant contributions at C(p^) 
(Eq. ([A16|) illustrates clearly the point). These redundancies also appear at higher orders if 



higher order tree Lagrangian terms are also included in the potential, V. 



B. OfF-shell versus on-shell 

Let us consider the BSE in the case of tttt scattering. The isospin amplitudes satisfy the 
symmetry properties given in Eqs. (|^) and (^^- Thus, they admit the following expansion 

TUp, fc) = E E E [^]^"^ ■ ■ ■ k'-^f' ■ ■ ■f''^ (36) 

where Tl^,„^^^.ur...vt,^ \P\ = Tyi...vi,^-tJiu..mi [P] and Ni and N2 run only over even (odd) natural 
numbers for the J = 0, 2 (/ = 1) isospin channel. In short hand notation we may write 

n{p,k)= y: t^mmI^i^^^V'^^ (37) 

Inserting this ansatz in the BSE, for simplicity let us consider first only those contributions 
where the free meson propagators (A*^*')) are used, we get 



Due to parity the number of indices in (a) and in {(3) ought to be either both even or both 
odd. Thus, we are led to consider the integral 

^......2„[p] .^.j _0_^V^...^M..A(o)(^^)A(o)(g_) (39) 

If we contract the former expression with P^^ we get, using that 

2^-^=(Tm^-Tm^)' (40) 



A{o)(g+) A{o)(g_)^' 
the identity 

P^Al-'^'-m = -i/ (04^''^'^ ■ ■ ■?'^^"AW(g+) (41) 
Shifting the integration variable 

P,J'2l-'''[P] = -i/ (04^^°Hg) ((g' - m") + (^ + m^) - g ■ P)' 

X (g - P/2Y'' ■■■{q- P/2Y''^ (42) 
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which corresponds to quadratic and higher divergences. As we have said before, only the 
logarithmic divergences should be taken into account since higher order divergences can be 
absorbed as redefinition of the renormalized parameters of the higher order terms of the 
Lagrangian. Thus, we would get that, up to these non- logarithmic divergences, the integral 
is transverse, i.e. when contracted with some P'^ becomes zero. The transverse part can be 
calculated to give 



+ Permutations 



+ P.D. (43) 



where C„fc(s) is a function]^ of s and P. D. means power divergences. Let us study now the 
integral l2n{.s) and use that 

= im' - |) + + A(°)(g-)-0 (45) 

then we get 

hn{s) = {m^ - + i / ^A(°)(g)g^"-^ (46) 

Again, the second term corresponds to power divergences. Applying this formula n times 

we get 

hn{s) = {m^-'-rh{s) + V.D. (47) 

All this discussion means that under the integral sign we may set = m? — s/4: and P-q = 
up to non-logarithmic divergences, i.e. 

i / ^F(g2;P-g)AW(g+)AW(g_)=F(m2-^;0)/o(.)+P.D. (48) 



Up to now we have considered the bare (free) meson propagator A(o). If we have the full 
renormalized propagator 

A{q) = [q^ -m^ -U{q^)y^ (49) 
where n(g^) is the meson self-energy with the on-shell renormalization conditions 

n{m^) = , n'(m2) = (50) 



^^It can be shown that 



Cnk{s) = (2n + l)!! ^^"+^^^^^- ^^^^ 



then we have the Laurent expansion around the mass pole, = rn^, 



A(g) = ^^ + + (51) 

The non-pole terms of the expansion generate, when applied to the kernel of the BSE, power 
divergences. Therefore, Eq. (^) remains valid also when the full renormalized propagator 
A(g) is used. 

In conclusion, the off-shellness of the BSE kernel leads to power divergences which have to 
be renormalized by a suitable Lagrangian of counter-terms (for instance the bare /'s— Lagrangian 
pieces at order 0{p^)), leaving the resulting finite parts as undetermined parameters of the 
EFT . Thus, if we iterate a renormalized potential we should ignore the power divergences 
and therefore ignore the off-shell behavior within the BSE and then solve 



T^ip, k) = Vlip, k)+ij j^JUq, k)A{q^)A{q_)Vl{p, q) (52) 
with 




^ ^ q'-iP-qY/s 

Thus, given the exact two particle irreducible amplitude, potential , the solution of Eq. (p, 
in a given isospin channel, and that of Eq. (|52D are equivalent. This important result will 
allow us in the next subsection to find an exact solution of the BSE given an exact two 
particle irreducible amplitude, V . It also might justify the method used in Refs. p, 
where the off-shell behavior of the Lippmann-Schwinger equation is totally ignored. 



C. Partial wave expansion 

We now show how the BSE can be diagonalized within this on-shell scheme for s > 0. 
To this end we consider the following off-shell "partial wave" expansions for Tp{p,k) and 
Vlip,k) 

oo 

TUp, k) = J2 (2-^ + P-p;k\P- k)Pjicos9,,p) (54) 

j=o 

and a similar one for Vp{p, k). The "angle" 6k^p is given by 

-p.k + {P-p)iP-k)/s 
cos9k,p ■= 1 r (55) 

[(p ■ p) Vs - 2 [(p .k)ys- P] 2 

and it reduces to the scattering angle in the CM system for mesons on the mass shell. For 
s > 0, the Legendre polynomials satisfy a orthogonality relation of the form 



i / ^A(g+)A(g_)P,(cos^,,,)P,,(cos^,,,) = 6jj,^^^^^Io{s) (56) 
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which is the generahzation of the usual one. Plugging the partial wave expansions of Eq. (|54D 
for the amplitude and potential in Eq. (|52D we get0 

Tu{p\ P-p;k\P-k) = Vu{p\ P-p;k',P-k) 

+ Tjj{m^ - ^ , ; P ■ k)Io{s)Vjj{p\ P-p;m'-^,0) (57) 

To solve this equation, we set the variable p on-shell, and get for the half off-shell amplitude 

Tjjim^ - f , ; P ■ k)-' = -/o(s) + V^jK -^,0;k\P- k)-^ (58) 

Likewise, for the full on-shell amplitude we get 

Tjj{m^ - £, ; - |, 0)-^ := Tu{s)-' = -Us) + Vjj{s)-^ (59) 

where V/j(s) = Vi.j{rn? — s/4, ; rn? — s/4, 0). Choosing the subtraction point sq we have a 
finite on-shell amplitude 

Tij{s)-^ - Tjj{so)-' = -(/o(s) - /o(so)) + Vjj{s)-' - Vuiso)-' (60) 

A finite (renormalized) inverse half-off-shell can now be obtained from the finite (renormal- 
ized) on-shell amplitude 

Tu{m' -'-,0-e,P- k)-' = Tu{s)-' + Vuim' - ^, ; fc^ P • k)-' - VUs)-' 

= Tij{so)-' - {Io{s) - /o(so)) 

+ (Vjj{m' - |, ; P ■ k)-' - Vjj{so)-'^ (61) 



Note that once we have the half off shell amplitude we might, through Eq. (pTf ), get the full 
off-shell one. The idea of reconstructing the full off-shell amplitude from both the on- and 
half off-shell amplitudes was first suggested in Ref. where a dispersion relation inspired 
treatment of the BSE, with certain approximations, was undertaken. 



D. Off-shell unitarity 

The off-shell unitarity condition becomes particularly simple in terms of partial wave 
amplitudes. Taking into account Eq. (^, that in this equation and due to the on-shell 
delta functions can be replaced hj in the arguments of the T— matrices, and taking 
discontinuities in Eq. (^) one can see that in terms of partial waves, off-shell unitarity reads: 

Tjj{p\P-p- k\P-k)-Tl{k\P-k-p\P-p) 

= Tu{m^ _ £, ; P ■ A;) Disc [Io{s)]T;j{m' - ^, ; P ■ p) (62) 



We preserve the ordering in multiplying the expressions in a way that the generalization to the 
coupled channel case becomes evident. 
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A solution of the form in Eq. (|57D automatically satisfies^ the off-shell unitarity condition 
of Eq. (|62|) . This property is maintained even if an approximation to the exact potential in 
Eq. ( p7D is made. 



E. Approximations to the Potential 

Given the most general two particle irreducible renormalized amplitude, Vp{p,k), com- 
patible with the chiral symmetry for on-shell scattering, Eq. ( |60D provides an exact T— matrix 
for the scattering of two pions on the mass shell, for any / J— channel. Eq. (|60D can be rewrit- 
ten in the following form 

TiAs)-' + Io{s) - VUs)-' = Tuiso)-' + Uso) - VrAsoy' = -Cjj (64) 

where C/j should be a constant, independent of s. Thus we have 

Tij{s)-^ = -Us) - Cu + Vu{s)-^ (65) 

or defining Wij{s, Cjj)"^ = —Cjj + Vjj{s)~^ the above equation can be written 

Tjj{sr' = -Io{s) + Wu{sr' (66) 



The functions Io{s) and Wjj{s) or Vjj{s) account for the right and left hand cuts |3! 
respectively of the inverse amplitude T/j(s)~^ The function Wij{s), or equivalently Vij{s), 
is a meromorphic function in C — R^, is real through the unitarity cut and contains the 
essential dynamics of the process. Thus, analyticity considerations ( |3^) tell us that the 
scattering amplitude should be given by 



pij^ 



Tf/is) = -lois) + (67) 



Q 

where P^'^(s) and Q''^{s) are analytical functions of the complex variable s G (D — R^. A 
systematic approximation to the above formula can be obtained by a Pade approximant to 
the inverse potential, although not exactly in the form proposed in Ref. p|. 



^^To prove this statement it is advantageous to note that above the unitarity cut, where the 
potential V is hermitian, 



r;^(A:2, P-k;p^P-p) = Vu{p\P -p^k^P-k) 



+ Vuim^ -^,0;k^P- k)rAs)Th{m' - £, ; P • (63) 
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where P^'^{x) and Qj/{x) are polynomials^ of order n and k with real coefficients. Crossing 
symmetry relates some of the coefficients of the polynomials for different /J— channels, but 
can not determine all of them and most of these coefficients have to be fitted to the data 
or, if possible, be understood in terms of the underlying QCD dynamics. This type of Pade 
approach was already suggested in Ref. []T2[. Actually, the results of the off-shell scheme 



presented in Sect. |T| are just Pade approximants of the type [n, k] = [1, 1] in Eq. (pq). 

A different approach would be to propose an effective range expansion for the function 
Wij{s)^ namely 



Wu{s) = Y.^ui^'^^'T^' (69) 

n=0 

or a chiral expansion of the type 

Wu{s) = T.^f3\j\s/m^)[^^j (70) 

where the (3^ij{s/m'^) functions are made of polynomials and logarithms (chiral logs). As 
we explain below, the ability to accommodate resonances reduces the applicability of these 
approximations (Eqs. (^) and (|70D) in practice. 

As discussed above, the function Wij{s) is real in the elastic region of scattering. Obvi- 
ously, the zeros of Tij{s) and Wij{s) coincide in position and orderQ. For real s above the 
two particle threshold, zeros could appear for a phase shift Sjj{s) = rni, G N, in which 
case a resonance has already appeared at a lower energy. For the present discussion we can 
ignore these zeros. The ffist resonance would appear at ^/s = for which 5ij{m^ = tx/2 
and hence ReTj}^(m^) = 0, i.e., — ReJo(mij) -|- W]~j{jnR) = 0, but ReJo(s) > 0, so a neces- 
sary condition for the existence of the resonance is that Wjj{mj^) > 0. On the other hand, 
the sign of Wjj{s) is fixed between threshold and the zero at Sjj{s) = nn. It turns out that 
due to chiral symmetry the functions Wu and Wqq are always negative for tttt scattering so 
the existence of the a and p resonances can not be understood. 

The above argument overlooks the fact that the change in sign of Wjj may be also 
due to the appearance of a pole at s = Sp, which does not produce a pole in T/j, since 
Tij{sp) = — l/Jo(sp). The problem is that if that Wu has to diverge before we come to the 
resonance, how can an approximate expansion, of the type in Eqs. ( |69D and (0), for Wjj 
produce this pole, and if so how could the expansion be reliable at this energy region ? To 
get the proper perspective we show in Fig. |^ the inverse of the function Wu (circles in the 
middle plot) for tttt scattering extracted from experiment through Eq. (|66|) . The presence 
of a pole (zero in the inverse function) in Wu is evident since the p resonance exists! 

We suggest to use Eq. (|65|) rather than Eq. (|66|) with the inclusion of a further unknown 
parameter Cjj, because in contrast to the function Wjj, the potential Vjj can be approxi- 
mated by expansions of the type given in Eqs. (|69|) and (ffOl). It is clear that with the new 



^'^Note that the first non- vanishing coefficient in Qn corresponds to the power (s — Am?Y . 

In TTTT scattering, the location of one zero for each / and J is known approximately (see e.g. 
Ref. |33i) in the limit of small s and m and are called Adler zeros. 
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constant C/j, KeTfj can vanish without requiring a change of sign in the new potential 
V/j(s). Thus, V/j(s) can be kept small as to make the use of perturbation theory credible. 
The effect of this constant on the potential extracted from the data can be seen in Fig. |^ 
for several particular values. 

After this discussion, it seems reasonable to propose some kind of expansion (effective 
range, chiral,. . . ) for the potential Vjj{s) rather than for the function Wij{s). The role 
played by the renormalization constant C/j will have to be analyzed. Thus, in the next 
subsection we use a chiral expansion of the type 

Vij{s) = Y.^ V}r\s/m^)^^j (71) 

where the ^^"^V/j(s/m^) functions are made of polynomials and logarithms (chiral logs). In 
this way we will be able to determine the higher order terms of the ChPT Lagrangian not 
only from the threshold data, as it is usually done in the literature, but from a combined 
study of both the threshold and the low-lying resonance region. 
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FIG. 3. Experimental isovector p— wave phase shifts (top panel), inverse of the functions Vii(s) (middle 
panel) and Kii{s) (bottom panel) as a function of the total CM energy y/s. Phase shifts are taken from the 
experimental analysis of Ref. |2^] and V{l^{s) is determined through Eq. ( |65| ) from the data of the top panel 
and using three values of the constant Cn: (Vn = Wn), —0.1 (as suggested by the results of Table |) and 
—0.038 (as suggested by the formula [log(rn//x) — l]/87r^ with a scale ^ of the order of 1 GeV, see Eq. (116)), 



which are represented by circles, squares and triangles respectively. Finally, Kn is also determined from 



the data of the top panel through Eq. (12C) 
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F. Chiral expansion of the on-shell potential 



The chiral expansion of the tttt elastic scattering amplitude reads 

TiAs) = Tf]{s)lf + Tjy + Tf]{s)lf + . . . (72) 

where t)j (s), t)j (s) and t)j (s) can be obtained from Refs. [Q, [|16| and ||35| respec- 
tively. For the sake of completeness we give in the Appendix ^ the leading and next-to- 
leading orders. 

If we consider the similar expansion for the potential Vjj given in Eq. ([Til) and expand 
in power series of / the amplitude given in Eq. (B^) we get 



P f 

+ ^ ( Vff + 2(/o + Cu) V^f V}f + (Jo + Cur [V}??) ■ ■ ■ (73) 
Matching the expansions of Eqs. (|72D and (|7^) we find 



2 Tr{2) rT,{2) 

m Vj\/ = T}j' 



Vif = Tff - (Jo + CuWil^r = - Cu[Tr]f 
with rjj defined in Eq. (|B^). From the unitarity requirement 



(74) 



Imri?(.) = -^^1 - ^[TS)(.)]^ . > (75) 

and thus, we see that Vjf is real above the unitarity cut, as it should be. 

Thus, in this expansion we get the following formula for the inverse scattering amplitude 

1 



Trj\s)IP + Trj\s)IP - (lois) + Cu)[T}y{sW/P + 

1 

TrAs)/p + Trj{s)/p~ cu[T\y{s)fip + 



Notice that reproducing ChPT to some order means neglecting higher order terms in 
s/ p and m?/p, which is not the same as going to low energies. With the constant Cu 
we may be able to improve the low energy behavior of the amplitude. Obviously, the exact 
amplitude Tij{s) is independent on the value of the constant C/j, but the smallness of 
Vij{s) depends on C/j. Ideally, with an appropriated choice of C/j, Vij{s) would be, in a 
determined region of energies, as small as to make the use of perturbation theory credible 
and simultaneously fit the data. 

Any unitarization scheme which reproduces ChPT to some order, is necessarily generat- 
ing all higher orders. For instance, if we truncate the expansion at fourth order we would 
"predict" a sixth order 
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fi'J = 2(/o + Cu) (rf; - Cjj[Ti^^r) T?) + (/o + C,jf[T?]f (77) 

and so on. 

Finally, we would like to point out that within this on-shell scheme crossing symmetry 
is restored more efficiently than within the off-shell scheme exposed in Sect. |T|. Thus, 
for instance, if we truncate the expansion at fourth order and neglect terms of order 1//^, 
crossing symmetry is exactly restored at all orders in the expansion (s — 4m^), whereas in the 
off-shell scheme this is only true if the u— and t— unitarity corrections are Taylor expanded 
around 4m^ and only the leading terms in the expansions are kept. 



G. Form Factors 



The interest in determining the half-off shell amplitudes lies upon their usefulness in 
computing vertex functions. Let Tp{p, k) be the irreducible three-point function, connecting 
the two meson state to the corresponding current. The BSE for this vertex function. Fig. ^, 
is then 

F^\k) = Tt{k) + ^ E i / yfri^pl?' kUa,^{q^mq.)T^p\q) (78) 

^ cd y^^r 

and using the BSE we get alternatively 




P/2-k,b p/2.q,d 

FIG. 4. Diagrammatic representation of the BSE type equation used to compute vertex functions. It is 
also sketched the used kinematics and a, 6, c, d are isospin indices. 



In operator language we have F = V + VGqF = F + TGqV, with Go the two particle 
propagator. The discontinuity in F is then given by the discontinuities of the scattering 
amplitude T and the two particle propagator Gq. 

The off-shellness of the kernel of Eq. ([78|) leads to power divergences which have to 
be renormalized by appropriate Lagrangian counter-terms. The resulting finite parts are 
undetermined parameters of the EFT (for instance /g, at order O(p^), for the vector vertex). 
The situation is similar to that discussed in Subsect. [IV A| for the scattering amplitude, and 
thus if we iterate a renormalized irreducible three-point function, F, we should ignore the 
power divergences and therefore ignore the off-shell behavior of the kernel of Eq. (|7Sp and 
then solve 
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^ cd y^^r 



(80) 



with q given in Eq. (0). The above equation involves the half-off or on-shell amphtudes for 
off-shell (/c^ 7^ m?) or on-shell {k"^ = m?) vertex functions respectively. 



1. Vector form factor 
The vector form factor, Fv{s), is defined by 

^7r'^(P/2 + k) T,\P/2 - A;) ^ (u^i^u - d^^d 



0) = -2tFvis)eab3h 



(81) 



where s = P^, a, b are Cartesian isospin indices, u, d and 7^ are Dirac fields, with flavor 
"up" and "down" , and matrices respectively. For on-shell mesons, k ■ P = and then the 
vector current | {u'^nU — d'^^d^ is conserved as demanded by gauge invariance. The vector 
form factor is an isovector, and we calculate here its /z = component, thus the sum over 
isospin in Eq. ( pOD selects the 7 = 1 channel. Then, we have 



Fy{s)k^ = Vy{s)k^ + \Vy{s) J ^Ti,(g, fc)A(g+)A(g_)r 
with rv(s) related with the irreducible vector three-point vertex by means of 



(82) 



(83) 



Using a partial-wave expansion for Tp{q, k) (Eq. (|5^) ) and the orthogonality relation given 
in Eq. ( |5^ ) the integral in Eq. ( p2D can be performed and for on-shell pions we obtain 



Fv{s) = Tv{s) + Tn{s)Iois)Tv{s) 



(84) 



The above expression presents a logarithmic divergence which needs one subtraction to be 
renormalized, choosing the subtraction point sq we have a finite form-factor 

(85) 
(86) 



Fv{s) = Tv{s) + Tn{s) [lo{s) + Cv) Tv{s) 
Cv = -Io{so) + 



Replacing Io{s) = —Tii^(s) — Cn + V^ii^(s) we get 



Ms) 



Vuis) 

Notice that Watson's theorem reads 



l + {Cv-Cu)Vu{s))Tv{s) 



(87) 



Fv{s-ie) Tii(s-ie) 



(88) 
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and it is automatically satisfied. That is a very reassuring aspect of the BSE approach. 
Indeed, in the literature it is usual to use Watson's theorem as an input, 



and employ it to write a dispersion relation for the form-factor. This procedure leads to 
so called Omnes-Muskhelishvili representation of the form-factor, which requires the 



introduction of a polynomial of arbitrary degree. In our case, not only the phase of Fv{s) is 
fixed, in harmony with Watson's theorem, but also the modulus, and hence the polynomial, 
is fixed at the order under consideration. 

The normalization of the form factor requires -FV(O) = 1, which allows to express the 
renormalization constant (Cy — Cn) in terms of Tn, Fy and Vn at s = 0, and thus we get 

Fv{s) = Tn{s) Tv{s) j-^ - + ^ ] (89) 

[Vu{s) Vu{0) Tii(0)ry(0) J 

From our formula it is clear that for s > 4m^, where both Ty and Vu are real, ReF(s)~^ = 
when ReTii(s)~^ = 0, in agreement with the vector meson dominance hypothesis. In 
Eq. (p9|). Til is obtained from the solution of the BSE for a given "potential", Vii, which ad- 
mits a chiral expansion, as discussed in Subsect. |1V F| . The irreducible three-point function, 
rv(s) admits also a chiral expansion of the type 

M.) = i + ^ + --- (90) 

Note that at lowest order, 0(p°), rv'(O) = -FV(O) = 1. The next-to-leading function r^(s) 
can be obtained if one expands Eq. (|89D in powers of 1//^ and compare it to the order 0{p^) 
deduced by Gasser-Leutwyler in Ref. [|T6|. Thus, we get 

rvW = i + {(i-£i)rm + s^(4-i)} + (91) 

For s » 4m^, Vii(s) and Fy(s) might increase as a power of s, whereas elastic unitarity 
ensures that Tii(s) cannot grow faster than a constant. Indeed, if Vn actually diverges 
or remains at least constant in this limit, then Tii(s) behaves like l//o(s), it is to say, it 
decreases logarithmically. To guaranty that the elastic form factor goes to zeroQ for s — oo, 
we impose in Eq. (|89D the constraint 

Fy(0) = l^ii(0)/Tii(0) (92) 



and thus, we finally obtain 



= IlMIM (93) 

Vuis) 



with Fy(s) given in Eq. (^) and F^(0) determined by the relation of Eq. (p2|). Thus, the 
vector form factor, at this order, is determined by the vector-isovector tttt scattering plus a 
new low-energy parameter: Iq. 



^^This behavior is in agreement with the expected once subtracted dispersion relation for the form 
factor fH 
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2. Scalar form factor 



The scalar form factor, Fs{s), is defined by 



7r"(P/2 + k) n\P/2 - k) (uu + dd 







(94) 



The scalar form factor is defined through the isospin-zero scalar source, thus the sum over 
isospin in Eq. ( ^OD selects the 1 = channel. Following the same steps as in the case of the 
vector form factor we get, after having renormalized. 



F.(s) 



d, 



Toojs) 
Voois) 
Cs — Coo 



l + dsV,^{s))Ts{s) 



(95) 



where Cs is given by Eq. (|5BD, with the obvious replacements V ^ S and 11 — > 00. Watson's 
theorem is here again automatically satisfied by Eq. (|95|). A similar discussion as in the case 
of the vector form-factor leads to set the ds parameter to zero, since the form factor goes to 
zero for s ^ oo [jT8| . Here again we propose a chiral expansion of the type of Eq. (ROl) for 



the irreducible three-point function, Ts{s 



p 



+ 



(96) 



and expanding the form factor in powers of 1//^ and comparing it to the order 0{p'^ 
deduced by Gasser-Leutwyler in Ref. [1161, finally get 



Fs{s) 
Fs{0) 

rf(o) 

To 
Fs{0) 



T sis) Toojs) VoojO) 
TsiO)TooiO) Voois) 



1 



1 + 



1 + 



rf(o) 



P To 

' 1 m'^Coo 



+ 



167r' 



:iU 



1 + 167r2Coo) 



To = 2B 
2b\i 



m 



167r2/2 



(^3 



1. 
2' 



(97) 



where is a low energy constant which measures the vacuum expectation value of the scalar 
densities in the chiral limit []TB[. Note that in contrast to the vector case the normalization 
at zero momentum transfer of the scalar form-factor is unknown. 



H. Numerical results: vrvr-phase shifts and form-factors. 



The lowest order of our approach is obtained by approximating the potential, 
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Vu - ^ Vi^ (98) 



with Vif given in Eq. (ffil). Thus, at this order we have three undetermined parameters 
Coo, C*ii, C*2o- At this level of approximation the d— wave phase-shifts are zero, which is not 
completely unreasonable given their small size, compatible within experimental uncertainties 
with zero in a region up to 500 MeV. Note that there is a clear parallelism between the C's 
parameters here and the J(f (4m^) low energy parameters introduced in Sect. |IT1[ Thus, the 



C's parameters will be given in terms of the U parameters, as we found in Subsect. [Ill F| , 
although some constraints on the I' s would be imposed since = = 0. However, one 
should expect more realistic predictions for the phase-shifts in the off-shell case, since there 
were more freedom to describe the data (four versus three parameters). For the sake of 
shortness, we do not give here any numerical results for this lowest order of the proposed 
approximation and also because they do not differ much from those already presented in 
Subsect. |111F[ The interesting aspect is that already this lowest order approximation is able 
to describe successfully the experimental phase shifts for energies above the certified validity 
domain of ChPT. We come back to this point in Sect. 0. In Subsect. [IV 1] we also investigate 



the convergence of the expansion for the potential, V , proposed in this work. Thus we will 
compare results at leading (discussed above) and next-to-leading (discussed below) orders 
for the function Wjj defined in Eq. (|66D. 

Further improvement can be gained by considering the next-to-leading order correction 
to the potential, which is determined by the approximation: 



Vu « ^ Vjf + ^ Vi9 (99) 



with v/j'^'' given in Eq. ([7^). At next-to-leading, we have nine free parameters (C/j, with 
IJ = 00, 11, 20, 02, 22 and k, i = 1, ■ ■ ■ A) for tttt scattering in all isospin channels and J < 2. 
Besides we have two additional ones, {B and Jq) to describe the scalar and vector form factors. 
As we discussed in Eq. (^7|), in this context the C's parameters take into account partially 
the two loop contribution, and thus they could be calculated in terms of the two loop ChPT 
low energy parameters. That is similar to what we did for the /^(4m^) parameters in 



Sect, m or what we could have done above with the C's parameters at the lowest order of 
our approach^. However, we renounce to take to practice this program, because of the great 
theoretical uncertainties in the determination of the needed two loop ChPT contributions: 



the two loop contribution have been computed in two different frameworks: ChPT ||35 



and generalized ChPT [^, and only in the former one a complete quantitative estimate 



of the low energy parameters is given. Furthermore, this latter study lacks a proper error 



analysis which has been carried out in Ref. |^^. On top of that a resonance saturation 



assumption |Q has been relied upon. Therefore, we are led to extract, at least, the C's 
parameters from experimental data. For the I' s parameters, we could either fit all or some 
of them to data and fix the remainder to some reasonable values of the literature, as we did 



^"^In both cases the unknown parameters could be calculated in terms, of the one loop ChPT low 
energy parameters, Vs. 
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in Subsect. [Ill F| . Here, we will follow a hybrid procedure. The justification will be provided 
a posteriori, since the error bars in the Vs are reduced in some cases. 



1. Electromagnetic pion form factor. 



Data on tttt scattering are scarce and in most of cases the experimental analysis relies on 



some theoretical assumptions because of the absence of direct vrvr scattering experiments K 
However, there are direct and accurate measurements on the electromagnetic pion form 
factor in both the space and time-like regions [^]. At next-to- leading order, the 

vector form-factor depends on Cu, h — h, U and Iq (see Eq. (|93|)). As we did for both sets 
of parameters in Eq. (|29|) we take 



L = 4.4 ±0.3 



(100) 



as determined by the the scalar radius 
|§ and ^ 



of Refs 
cuts. 



I and fit the another three parameters to the data 
Results can be seen both in Fig. ^ and Table y, for different energy 



h = 4.4 


^ < 483 


^/s < 700 


^ < 808 


< 915 


^ < 1600 


h - h 
Cn 
k 

xVdof 
num. data 


-6+1 

-n 1 5+0-0'' 

^•-'-'3-0.05 

19.0 ±0.3 
1.2 

63 


-6 1+0-1 
-0.112 ±0.010 
19.14 ±0.19 
1.2 
76 


_r: 02+0-03 
'^•^^-0.04 

-0.1099 ±0.0018 
20.84 ±0.06 
3.9 
112 


-5.93 ±0.03 
-0.1098 ±0.0018 
20.82 ± 0.06 
3.5 
126 


-5.861 ±0.010 
n iiqt^+0 0016 

-U.iiy£l_Qooi5 

21.09 ±0.06 
3.3 

184 



TABLE II. Results of best fits of the next-to-leading approach (Eqs. (p^), (93) and (99)) to the electro- 
magnetic pion form factor. Data are taken from Refs. (space-like, — (503MeV)^ < s < — (122.47MeV)^ 
) and (time-like, 320 MeV < ^/s < 1600 MeV). Errors in the fitted parameters are statistical and have 
been obtained by increasing the value of by one unit. We fix I4 = 4.4 and show the variation of the fitted 
parameters, their statistical errors and x^/dof with the used energy cut (given in MeV) for the time-like 
region. We also give, in the last row the number of piece of data of each of the best fits. 
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h-h 


Cn 


k 


4.1 


-6 2+0-1 


-0 112 + 0011 
'-'•1-1^ - 0.010 


19.13 ±0.19 


4.4 


-6 1+0-1 


-0.112 ±0.010 


19.14 ±0.19 


4.7 


-6 1+0-1 

"•1-0.3 


-0.111 ±0.010 


19 15 + 0-19 



TABLE III. Results of best fits of the next-to- leading approach (Eqs. (p5|), ( p3[ ) and (p9|)) to the 
electromagnetic pion form factor for different values of 1^. Data are taken from Refs. and |4^. In all 
cases, the fit range is —503 MeV < s/|s|5 < 700 MeV. Errors in the fitted parameters are statistical and 
have been obtained by increasing the value of by one unit. 
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FIG. 5. Best fits of tlie next-to- leading approach (Eqs. (|65| ) and (p3|)) to the electromagnetic pion form 
factor. Data are taken from Refs. [|5| (space-like) and |Q (time-like). Results obtained with the parameter 
set determined by the entry ^/s < 700 {y/s < 1600) of Table || are displayed in the left (right) panels. 

As can be seen both in the last column of Table y and in the right plots of Fig. a 
fairly good description of the data from —500 MeV < s/|s|2 < 1600 MeV can be achieved. 
However, our aim is not only to describe the data but also to determine the low energy 
parameters which determine the chiral expansion of both the electromagnetic form factor 
and the vrvr-amplitude. If one look at the values of /dof quoted in Table |I|, one realizes 
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that there is a significant change (1.2 versus 3.9) when data above 700 MeV are considered. 
Besides, the fitted parameters turn out to be not statistically compatible below and above 
this cut in energies. Although with a small variation on li — I2, and Cu one still gets a 
reasonable description of the form-factor in the whole range of s, one can appreciate, in the 
two last rows of plots of Fig. |^, that the higher the energy cut, the worse the description 
of the low energy region (—500 MeV < s/|s|2 < 600 MeV). Indeed, this is the region 
where the chiral expansion is supposedly more trustable. Neither higher orders in the chiral 
expansion proposed in this work, nor the effect of non-elastic channels (like 47r and KK 
production) have been included at this level of approximation, and both might account for 
the discrepancies appreciated in the top left panel of Fig. ^ at high energies. The small 
changes in the fitted parameters above and below the energy cut of 700 MeV effectively 
incorporate these new effects into the model. Thus, our best determination of the one-loop 
ill — I2 and Iq) and two-loop Cu low energy parameters come from the best fit with an 
energy cut of 700 MeV in the time-like regionQ. Besides the statistical errors quoted in 
Table ||, we should incorporate some systematics due to the uncertainties in 1^. From the 
results reported in Table |T| we deduce that such systematic errors are much smaller than 
the statistical ones quoted in Table || and can be safely ignored^. Thus, we get 

li-l2 = -Q.ltol, Is = 19.14 ±0.19, = -0.112 ±0.010 (101) 

which provides us with an extraordinarily precise determination of the difference li — I2 and 
of Zg- The authors of Ref. [T^ have completed a two loop calculation for the scalar and 
vector form factors, finding Iq = 16.0 ± 0.5 ± 0.7 which differs by as much as two standard 
deviations from our result. 



2. Elastic TTir— scattering in the p— channel. 

The parameters of Eqs. (|100| ) and ( |101|) uniquely determine, at next-to-leading order in 
our expansion (Eq. (|99D ), the vector-isovector vrvr scattering phase shifts. Results are shown 
in Fig. ^. The solid line has been obtained with the central values quoted in Eqs. (|100|) 
and (|101|) . There also we show the 68% confidence limits (dashed lines), obtained by as- 
suming uncorrelated Gaussian distributed errors in the parameters quoted in Eqs. ( |100| ) 
and (|101|) and propagate those to the scattering phase shifts by means of a Monte Carlo 
simulation. For the difference /i — I2 we have assumed a symmetric error of magnitude 0.3. 



^^Lower cuts lead to similar x^/dof and statistically compatible fitted parameters, but with 
appreciably higher errors, as expected from the reduction on the amount of experimental data 
points. Obviously, we take the highest cut-off, 700 MeV, which keeps the \^ /dof around one, 
maximizes the number of data and hence minimizes the errors on the fitted parameters. 

^^Vector form factor data are almost insensitive to /4 which has prevented us to determine it from 
a four parameter fit, but rather to fix it to the value obtained in Ref. p8| . 
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FIG. 6. On-shell next-to-leading order BSE (Eqs. (|6^) and (|99[)) vector-isovector phase shifts as a func- 
tion of the total CM energy y/s. The solid line has been obtained using the central values quoted in 



Eqs. (100) and (101). Dashed lines are the 68% confidence limits assuming Gaussian and uncorrelated error 
distributions for the parameters. For the difference li — I2 we have assumed a symmetric error of magnitude 
0.3. Circles stand for the experimental analysis of Refs. ||2^] and and the triangles stand for the Frogatt 
and Petersen phase-shifts (Ref. ||2^) with no errors due to the lack of error estimates in the original analysis. 



Besides we also findQ 



rrir 



76Atll [MeV], Tp = U9tl^ [MeV] 



(102) 



in fairly good agreement with experiment and we also obtain the threshold parameters 
reported in Table |V|. 



3. Elastic TTTT— scattering for s— and d— waves. 



In the original work of Gasser-Leutwyler [|r^ and in the subsequent analysis of Ref. 
(set II in this latter reference) the parameters li^2 are measured through the c?— wave scat- 
tering lengths. Both scattering lengths and specially the isotensor one, suffer from large 
experimental uncertainties [^] which lead to quite big errors on li^2- Here, we have in 
Eq. ( |101|) a precise determination of the difference h — h, thus we can combine it with the 
experimental measurement of to determine both parameters /i 2- In this way we avoid 
to use a22, whose relative error is about 13 times bigger than that of ao2- We get: 



/i = -0.5 ±0.5, /a = 5.6 ±0.5 



(103) 
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We use the procedure described in Subsect. IIIF. 
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where we have used ao2 = (17±3) -10 [0]. The error on the scattering length dominates 
by far the errors on li 2- Fixing from the SU{3) mass formulae 

h = 2.9 ± 2.4 (104) 

to complete the analysis of the tttt— scattering at next-to-leading order in our expansion 
(Eq. (pUD), we still have four undetermined parameters ( Cjj, with IJ = 00, 20, 02, and 22). 
We determine these four parameters from four independent best fits to the data, fixing the 
li parameters to the central values given in Eqs. (|TUD|),(|TU3|) and ( p.04| ). Results are reported 
in Table There, we give Cjj and their statistical errors which have been obtained by 
increasing the value of the corresponding by one unit. In addition to these statistical 
errors, we have some systematic errors due to uncertainties in the one loop parameters Tj. 
One way of taking into account both types of errors would be to perform a simultaneous 
eight-parameter (C/j, with IJ = 00, 20, 02, and 22 and /i.2,3,4) fit to data. However and 
due to the limited quality of the available experimental data and the multidimensional 
character of the fit, it is very difficult to single out a stable minimum advising against 
carrying out this procedure directly. Thus, we have designed an alternative strategy to 
estimate the systematic uncertainties in the C's parameters induced by the errors in the 
one loop parameters. We fix central and statistical errors of the parameters /i, 2,3,4 to the 
values quoted in Eqs. ( p.O0|) ,( p!O^) and (|104D , and generate a sufficiently large sample of low 



energy parameters ^1,2,3,4 randomly distributed according to a Gaussian. For each set of four 
Vs parameters we fit, in each angular momentum and isospin channel, the Cjj parameters 
to the experimental phase shifts. In this way we generate distributions for each of the C's 
parameters. For the central values of the C's parameters we take the results of the fits 
obtained with the central values of Eqs. ( |100| ),( p!03| ) and ( |104| )p^. The 68% confidence limits 
of each of the distributions give an estimate of the systematic errors in the C's parameters. 
Since the out-coming C's parameter distributions are, in general, not Gaussian we end up 
with asymmetric errors. Results of this procedure are also reported in Table 

A final detail concerns the choice 610 MeV as the energy cut in the scalar-isoscalar 
channel. This is justified to avoid any possible contamination from the KK channel in the 
sub-threshold region or from higher chiral orders, not included in neither case, at this level 
of approximation. 

Solid lines in Fig. |^ are the predictions for the phase shifts in each angular momentum- 
isospin channel. At each fixed CM energy, our prediction for the phase shift will suffer 
from uncertainties due to both the statistical errors on k and Cjj parameters and also the 
systematic errors on the latter ones induced by the statistical fluctuations of the former ones. 



This is not the same as computing the mean of the generated distribution. Indeed, the mean 
values in general can not be obtained from any specific choice of the input parameters, in this case 
the Vs. This is the reason to quote central values from results obtained from central values of the 
input parameters. Furthermore, the statistical errors obtained before by changing the by one 
unit, are referred to the values resulting from best fits to the data with input parameters fixed to 
their central values. Both choices do not differ much as long as we are dealing with relatively small 
asymmetries in the distributions. We should stress that, what it is statistically significant is the 
68% error band rather than the choice of the central value. 
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as we discussed above. We use a Monte Carlo simulation to propagate independently both 
sources of errors and finally we add in quadratures, respecting the possible asymmetries, 
the statistical and systematic phase shift errors. The dashed lines in Fig. |^ join the central 
value for the phase-shift plus its upper total error or minus its lower total error for every 
CM energy. We see that, in general, the agreement with experiment is good except for the 
scalar-isoscalar channel, where one clearly sees that there is room for contributions due to 
the KK inelastic channel above the energy cut (610 MeV) used for the best fit. 



IJ 


00 


20 


02 


22 


Cij 
X^/num. data 


n 099+0-001 +0-005 
^•'-'^^-0.001-0.003 

34.7 / 20 


-0.058 ± o.oo2l:[}:°25 

20.5 / 21 


-0.203 ± 0.003;[]:|j^[] 
11.4 / 24 


-0.5±0.2l°;| 
0.4 / 24 



TABLE IV. Next-to-leading (Eq. (|9|)) parameters fitted to the experimental data of Refs. ||2| - (26) 
(with an energy cut of 610 MeV) for / = J = 0, of Refs. |3 and @ for / = 2, ./ = and finally of Ref. (2| 
(with an energy cut of 970 MeV) for the d— wave channels. In the latter case and due to the lack of error 
estimates in the original references, we have assumed an error of 0.4 in all phase shifts, this assumes that 
the errors of the d— wave phase-shifts given in Ref. |23| , affect only to the last digit. The central values and 



statistical errors of the parameters /i,2,3.4 have been fixed to the values quoted in Eqs. ( |100| ),(103) and (104). 
We give two sets of errors in the fitted parameters. The first set corresponds to statistical errors and have 
been obtained by increasing the value of the corresponding 'x^ by one unit. The second set are systematic, 
induced by the uncertainties in ^1,2,3.4, and have been estimated by means of a Monte Carlo simulation, see 
text for details. 



IJ 


00 


11 


20 


02 


22 




216+°-°°^ 


0.0361 ± 0.0003 


-0.0418 ± 0.0013 


input 


0.00028 ± 0.00012 


(exp) 


0.26 ±0.05 


0.038 ± 0.002 


-0.028 ±0.012 


0.0017 ± 0.0003 


0.00013 ± 0.00030 




0.284i°;°?l 


0.0063 ± 0.0005 


-0 07^+0-002 


481 
201600 TT^/-* 


277 
201600 7r^/4 


(exp) 


0.25 ±0.03 




-0.082 ± 0.008 







TABLE V. Threshold parameters a/j and 6/j (obtained from an expansion of the scattering ampli- 
tude j2^, ReT/j = — 167rm(s/4 — inn?')'-'[ai,j + bjj{s/4 — w?) + • - -] close to threshold) deduced from the 
0{p'^) results reported in Table. |^ and Eq. ( |101| ). For all channels, except for the I — J — 1 one, the errors 
have been obtained by adding in quadratures those induced by the uncertainties in the /'s parameters and 
those induced by the experimental errors in the tttt phase-shifts. For the vector-isovector channel, the error 
treatment is simpler thanks to the simultaneous determination of — ^2 and Cn from a best fit to data 
and to the negligible effect of the error bars of Z4 on the latter fit. Thus in this channel, uncertainties have 
been obtained by assuming uncorrelated Gaussian distributed errors in the parameters quoted in Eqs. (IOC) 
and (101). We also give the known experimental values compiled in Ref. |^9t] . 
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FIG. 7. Several tttt phase shifts as a function of the total CM energy y/s obtained within our 
next-to-leading approximation (see Eq. (|99|)). The top (bottom) panels correspond to / = (2) whereas the 
left (right) panels correspond to J = (2). Experimental data are taken from Refs. ~ |^ (/ = J = 0) 
from Refs. |2^] and [Q {I — 2, J — 0) and finally from Ref. |Q for the d— wave channels. Triangles stand 
for the Frogatt and Petersen phase-shifts (Ref. 123]) with no errors due to the lack of error estimates in 



the original analysis. Solid lines have been obtained with the central values quoted in Table IV and in 



Eqs. (|l00| ),(103) and (104). Dashed lines show our total uncertainties, statistical and systematic, added in 
quadratures, see text for more details. 

In Table we also give the deduced threshold parameters. We propagate both types 
of errors discussed above to these observables, by means of a Monte Carlo simulation, and 
finally we add both errors in quadratures. In general, for those parameters which have been 
measured, we find agreement within experimental uncertainties, though our final errors are 
always significantly smaller than the experimental ones compiled in Ref. ||49|| . 
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4- Numerical comparison with two loop ChPT. 



Set 


(tree) +(lloop) +(21oop) 


total 


experiment 


I 

TT 

aoo il 

BSE 


0.044 ± 0.005 0.016 ± 0.003 
0.156 0.039 ± 0.008 0.013 ± 0.003 

0.045 ± 0.004 0.011 ; 


0.216 ±0.009 
0.208 ± 0.011 

216+°-°°^ 

U.^iU_Q QQg 


0.26 ± 0.05 


I 

U TT 

BSE 


0.069 ±0.010 0.027 ±0.007 
0.179 0.059 ± 0.024 0.019 ±0.011 
0.072 ± 0.009 0.006^[j:[![]^ 


0.275 ±0.016 
0.256 ± 0.034 

0.284i°:°?l 


r\ cr 1 r\ no 

0.25 ± 0.03 


I 

1 r\ TT 

10 • ail ■I-'- 

BSE 


0.073 ±0.010 0.025 ±0.006 
0.297 0.058 ± 0.033 0.018 ± 0.005 
0.064 ± 0.003 


0.395 ±0.014 
0.374 ± 0.034 
0.361 ± 0.003 


n 00 1 n no 

0.38 ± 0.02 


I 

1 r> I. TT 

10 • Oil ■I-'- 

BSE 


0.048 ± 0.006 0.031 t 
0.034 ± 0.033 0.020 I QQQg 
0.0389 ±0.0017 0.021 ±0.003 


0.080 _ 0.009 
r\ n c A 1 n nor\ 

0.054 ± 0.029 
0.063 ± 0.005 




I 

1 r\ TT 

10 • a2o II 

BSE 


0.028 ±0.018 0.004 ±0.002 
—0.446 0.008 ± 0.031 0.000 I qqq3 
0.028 ± 0.016 0.000 t 


-0.414 ± 0.020 
— 0.438 ± 0.032 
-0.418 ±0.013 


n 00 1 n 1 

—0.28 ± 0.12 


I 

^0 ■ hnn TT 

BSE 


0.17 ±0.04 0.01 ±0.01 
—0 89 1 ± 05 00 ± 01 
0.16 ±0.04 O.Oll[];[!^ 


-0.72 ± 0.04 

-0 75+0-02 


_n 89 ± 08 


I 

102 . ao2 II 

BSE 


0.181 ±0.025 0.079 ±0.016 
0.117 ±0.026 0.053 ±0.018 
0.170 ±0.030 


0.260 ± 0.036 
0.170 ± 0.030 
input 


0.17 ±0.03 


I 

10^ • 022 II 

BSE 


0.21 ±0.13 -QM^QM 
0.12 ±0.44 O.Oll[];J^ 
0.28 ±0.12 


0.20 ±0.10 
0.13 ±0.30 
0.28 ±0.12 


0.13 ±0.30 



TABLE VI. Separation in powers of 1//^ of some of the threshold parameters reported in Table 0. 
For comparison we have also compiled the results reported in Table 1 of Ref. which are obtained 
from the two-loop analysis of Ref. [^5| supplemented with proper error estimates. The Sets I and II 
refer to those define in Ref. The authors of Ref. assume some resonance saturation to give 

numerical values to the 0{p^) parameters at certain scale around 750 MeV (see that reference for details). 
In addition, Set I uses: = -1.7 ± 1.0 , h = 6.1 ± 0.5 , h = 2.9 ± 2.4 , [4 = 4.3 ± 0.9 and Set II uses: 
h = -0.8 ± 4.8, ^2 = 4.45 ±1.1 ,[3 = 2.9 ±2.4, [4 = 4.3 ±0.9. On the other hand, BSE stands for the results 
of the present work at next-to- leading order, with parameters given in Table ^ and Eqs. (IOC), (101), (103), 
and (104). Tree, one-loop, two-loops and total labels stand for the threshold parameters calculated at 0{p^), 
0{p'^), 0{p^) and all orders, respectively. We also give the known experimental values compiled in Ref. 
All threshold parameters are given in pion mass units. Note, that due to correlations between the several 
orders, in the 1//^ expansion, contributing to the threshold observables, the errors of the total quantities 
can not be simply added in quadratures; we use a Monte Carlo simulation to propagate errors. 
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A very interesting feature of the present treatment at next-to-leading order is that having 
fitted the C's to the data, we can use those parameters to learn about the two and even 
higher loop contributions (see Eq. (|T^)). The separation of the amplitudes near threshold, 
in powers of 1//^, is presented in Table ^ and compared to recent two loop ChPT calcula- 
tions supplemented with proper error estimates [^]. Despite of the resonance saturation 
hypothesis[3 assumed in Refs. [^,0] , in general we see that the BSE predictions are not 
less accurate than those of these references. Both sets of results, BSE and ChPT, suffer 
from systematic errors induced by the higher order contributions. Those are not included 
in either case. 



5. Scalar pion form factor. 



To end this subsection, in Fig. ^, we present results for the scalar form-factor (Eq. 
normalized to one at s = and computed at next-to-leading order accuracy, within our 
approach. It is to say, it uses, Voo(s) given in Eq. (|99|) with the parameters presented in 
Eqs. ( |10(J| ),( [TUB| ) and ( |104| ) and Table |V[ We also compare our results with those obtained 



in Ref . . Our results disagree with those of Ref . above 500 MeV, presumably due to 



the role played by the sub-threshold KK effects. 



I. Comparison of leading and next-to-leading on shell BSE predictions. 

In this subsection we compare results for elastic tttt— scattering at leading and next-to- 
leading accuracy within the on-shell BSE scheme. This kind of comparison is not generally 
undertaken in the literature regarding unitarization methods, where the main goal is just to 
fit the data, rather than study also the convergence of the expansion. Since in our framework, 
the potential plays the central role, such a comparison is naturally done in terms of it or in 
terms of the function Wjj{s) = Tfj{s) +/o(s) defined in Eq. (|5^. These Wjj functions can 
be directly extracted from dataQ in order words, they can be derived from experimental 
phase-shifts. In Fig. ^ we present experimental, BSE-leading and -next-to-leading results^, 
together with both experimental and theoretical error estimates, for Wqq^, W2o^, W02, W22 



Such an approximation reduces a priori the errors on the two loop contributions with respect 
to the present framework where these higher order corrections have to be fitted directly to experi- 
mental data. 

^^Note, the functions Vij can not directly be obtained from experiment, because of the unknown 
constants Cjj. 

^^Note, that because the presence of resonances in the a— and /)— channels, the corresponding 
ly— functions should have a pole (see discussion related to Fig. so we present the inverse for 
these channels and also, for a better comparison, for the isotensor-scalar one. On the other hand 
for d— waves the lowest BSE order approximation gives T = 0, and therefore for those channels 
we present W itself. 
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FIG. 8. On-shell next-to-leading order BSE (Eqs. l pT\ ) and (|99|)) scalar form factor (circles) as a function 
of s. The used parameters are those given in Eqs. (|10|),(|10|) and ([l04|) and Table 0. For comparison we 



also show the results of Rcfs. |40 (solid line) 



and W^l . Next-to- leading results are obtained from Eq. (|99D . Central values and uncer- 
tainties of the needed low energy constants were discussed in Subsect. p.V H| (Table ^ and 
Eqs. ( |100|) , (|101|) , ( |103| ), and ( p.04| )). Errors have been propagated via a MonteCarlo simula- 
tion and statistical and systematic errors have been added in quadratures. Leading results 
have been obtained from Eq. (|98D where the unknown low energy constants Cqo, Ch and 
C20 have been obtained from three one parameter best fits to the experimental data, since 
at lowest order the three channels are independent. For the s— wave channels we have fitted 
the same set of data and use the same energy cuts as in Table ^ and Fig. |^. We have 
obtained 



00 



-0.0273 ± 0.0004 xV^um. data = 31.2/20 



(105) 



C20 = -0.051 ± 0.014 xVnum. data = 21.4/21 

For these two channels both leading and next-to-leading approximations provide similar 
descriptions of the data. For the p— channel the situation is different. We have fitted Cu to 
the data of Ref. ||2T[. Fits from threshold up to 0.9 GeV give values of x^/dof of the order 
of 12, being then highly unlikely that such a big discrepancy between theory and data is 
due to statistical fluctuations, disqualifying any error analysis. We find significantly lower 
values of x^/dof for smaller energy cuts, being the optimum choice obtained with a energy 
cut of about 690 MeV, for which we find: 

Cn = -0.1275 ±0.0006 x^/dof = 4:.l (106) 

The above error is quite small, due to large value of x^- Leading and next-to-leading esti- 
mates for all three parameters Coo, C20 and Cu are compatible within two sigmas. 
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FIG. 9. Experimental (circles), BSE-leading (dashed lines) and -next-to-leading (solid lines) results, to- 
gether with both experimental and theoretical (dash-dotted and dotted lines for leading and next-to-leading 
predictions respectively) error estimates, for Wg^""^, W^^, W02, W22 and W^-^ as a function of the CM energy. 
In some cases the theoretical errors can not be seen due to their smallness. See Subsect. IV I for more details. 



41 



Results from Fig. ^ show a good rate of convergence of the expansion proposed in this 
work. 



V. THE NON-PERTURBATIVE NATURE OF THE C/j PARAMETERS. 

In this section we hnk the renormahzation Cjj parameters to some non-perturbative 
effects of the underlying theory, such as the the existence of resonances. 

For simplicity, let us consider the p— channel (/ = J = 1) for which the lowest order 
inverse amplitude, within the on-shell scheme (Eq. (P5|)), reads 

(107) 



where we have chosen the subtraction point Sq = — //^ . From above, we identify Cn by 

Op 
fi'^ + Am 

which should be independent of the scale /i^. Thus we have 



Cu = -Tn^( V) - /o(-/i^) + (10^ 



T-(.) = -(/o(3) + C^n)-^3^ (109) 
At the resonance, s = m^, ReTii(m^)~^ = 0, and thus if we look at the real parts, we have 

Cn = -Re/oK^ ^ 



m'i — Aw? 



1 ,og!i-54 + o(!!4) (110) 



-0.110 



in qualitative good agreement with the findings of previous sections (see, for instance, the 
Iq entries for this channel in Table | or Eq. ( |101j ) ) . 
Similarly for the scalar-isoscalar channel, 

T„-„'(.) - T-( V) = - (/o(.) - /o(-.^)) - P (j^^ + (111) 

where we have chosen the subtraction point Sq = — /i^ and it can be different for each IJ 
channel since we may take Tij{—ix^) as the fitting and a priori unknown parameter. As 
before, we can identify Coo 

Cm = -Too^( V) - /o( V) + (112) 

From all our previous discussions and fits it is rather clear that in the / = J = channel the 
a resonance cannot be completely understood without inclusion of sub-threshold KK contri- 
butions. Nevertheless, let us define the point s = m„ as the one fulfilling ReToo(m^)~^ = 0, 
yielding 
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' ,ogi-4 + 0(^) (113) 



f« -0.043 

for m„ = 600 MeV. This is not completely unreasonable since the phase shift at about 
s = 600MeV go up to 60°, which is equivalent to take ReToo(m^)~^ small. From this 
discussion, we understand the origin of the significant difference of size between Cu and 
Coo- Indeed, neglecting the logarithmic contributions, 

which clearly points out that the reason for the different size between the two C parameters 
relies on the fact that the potential in the p— channel is about six times smaller than that in 
the scalar-isoscalar one. This has also been pointed out by the authors of Ref. [0] who also 
use this argument to stress the dominant role played by unitarization in the a channel. The 
authors of this same work, also claim that the quite different values of Coo and Cu could 
not be understood if the same scale or cutoff, as they used in a previous work 0, are used 
for both channels. Indeed, if one uses dimensional regularization to compute the divergent 
integral Io{s) and in the MS scheme, one gets the finite piece 

1 f 2 J ^ 

167?^ 1 ^1/2 



^^(^) = T7-I -2 + log- +/o(.) (115) 



being p some renormalization scale. Thus, one might identify 

for all angular momentum and isospin channels. The above equation gives values for the C 
parameter of about —0.04 when scales, u, of the order 1 GeV are used. This agrees reasonably 
well with the values found for Coo, but to get values of about —0.11, as it is the case for Cu, 
one needs unrealistic scales or cut-offs of the order of 300 GeV. All confusion, comes from 
the interpretation of the C parameter as a renormalization of the bubble divergent integral 
/o(s), instead of a renormalization of the full amplitude once that a renormalized potential 
is iterated by means of the BSE. A direct, comparison of Eqs. ( |116D and ( |10(j| ) or ( |112[ ) 



shows that the free parameters Tjjl—fi^) allow for the use of the same value of the scale n 
to generate such different values for the C/j parameters in the scalar-isoscalar and vector- 
isovector channels. Furthermore, as we said above, the C parameters defined in Eq. (|108|) 
or (|112|) do not depend on scale. On the light of this discussion, the numerical coincidence 



between the values provided by Eq. ( |116| ) and Coo might be considered accidental. 

From the former expressions and neglecting corrections of the order of the (9(m^/m^^) 
we find simple expressions for the amplitudes in the two channels considered in this section, 

Tn\s) ^ --i^log(-4) - 6/'(- - 
167r^ rrip s 

Too\s) ^ -T^log(-4) - f{- - ^) (117) 
167r^ mi s mi 
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where we take log(— s/m|.) = log(|s|/m|j) — i7r6(s). The above expressions lead to approx- 
imate expressions for the decay widths (see Eq. (|30D ) 



_ ml 1 



967r2/2 

1 

IGn'^P 1 - 



(119) 



167r2/2 



For the p— channel, this corresponds to a coupling (^ptttt ~ (^p/v^/)) the value predicted 
by the KSFR relation 0. KSFR predicts, for = 770 MeV a width of Tp = 141MeV, in 
excellent agreement with the experimental value. Such a derivation of the KSFR relation is 
actually not completely unexpected since this relation is a consequence of PCAC plus vector 
meson dominance through current-field identities; in our model vector meson dominance is 
realized, i.e. the real parts of the inverse scattering amplitude in the I = J = 1 channel and 
of the inverse vector form factor vanish simultaneously, as we pointed out after Eq. (pQ]). 
The width for the a meson turns out to be very large making this way of determining it 
doubtful. Notice that in both cases the conditions for non negative widths are < QGvr^/^ 
and < IGvr^/^ respectively. 



VI. COMPARISON WITH OTHER UNITARY APPROACHES 

As we have hopefully pointed out along this work, the unitarization generated through a 
solution of the BSE becomes extremely important to describe the resonance region and also 
"guides" the theory in the low energy regime, allowing for a reasonably accurate description 
of low energy and threshold parameters. 

In what follows we want to compare the BSE with some other approaches pursued in 
the literature for on-shell amplitudes which also incorporate elastic unitarity. Some of the 



methods presented here have already been discussed, |5^, in the chiral limit. 

All unitarization schemes provide the same imaginary part for the inverse amplitude and 
they differ in the way the real part is approximated. Thus, it is not surprising that formal 
relations among them can be found. 

A. K-matrix 

A usual unitarization method consists in working with the 7^— matrix, defined as 
KYj = Re{Tr/)=Tr/-zp 



p = -ImJo(. + ie) = \' ' ^ = —\1 (120 

In the bottom graph of Fig. |^, we plot the inverse of the i^'— matrix for the isovector p— wave 
channel. The presence of a pole (zero in the inverse function) in Ku limits the use of any 
"perturbative" expansion of Ku in the region of the resonances and this fact prevents the 
use of this unitarization scheme above threshold when there exists a resonance. 
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B. Dispersive Representation of the Inverse of the Partial Wave Amplitude. 



Partial wave amplitudes, T/j, have both right and left hand cuts. Besides, the only poles 
that are allowed on the physical s sheet are the bound states poles, which occur on the 
positive real axis below 4m^. Poles on the real axis above threshold violate unitarity, and 



those on the physical sheet that are off the real axis violate the Mandelstam hypothesis ||32 
of maximum analyticity. Thus, Tjj and Tfj have similar right and left cuts. However, 
partial wave amplitudes have zeros, which in turn means that their inverses have poles on 
the physical sheets. 

A seemingly significant advantage of considering the dispersive representation of the 
inverse of a partial wave amplitude is that its right hand cut discontinuity is just given by 
phase space in the elastic region, thanks to unitarity. The relevance of this property was 
already emphasized in Ref. [^]. For simplicity, let us suppose that Tjj, has a zero of order 
1 for s = Sj^. 



If we define, 



r= \im{s-SA)Trl{s) (121) 



S^SA 



assuming that \Tfj(s)\ < s as s ^ oo then each inverse amplitude satisfies a once-subtracted 
dispersion relation p, with s = si, the subtraction point (dropping the labels I, J for 
simplicity) 

\{s-Sa) {Si-Sa)J TT J'km^ [S' - Si){s' - S) 

{S — Si) ds' ^ I 



Elastic unitarity (Eq. (|26|)) allows to perform the integral over the right hand cut, thus one 
gets 

T'\s) - T-\s,) = (-^ - r^-^] - {m - Us,: 

\{s-sa) {si-sa)J ^ 

is-Sl) ds' ^ rr.-1/l- 



+ ^I., (.--..)(^'-^) ""^ ' 

Then one clearly sees that all the dynamics of the process is contained in both the poles of 
T^^ (zeros of T) and in the left hand cut. This latter contribution can be split as follows: 

[s-si) /-o ds' , 1 



ImT-^(s') = — — (124) 



vr J-oo{s' - si){s' - s) ' ' V(s) V(si) 



^For the s— waves, these are demanded by the Adler condition, for the p— waves, the zero is 
kinematic being at threshold, sa = 4m^. In both latter cases, the zeros are of order one, higher 
partial waves have higher order zeros. 
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where 



where P is a renormahzation constant, though in general might be a polynomiaL Thus, 

T^^{s) can be written as 

T-\s) = -Io{s)-k + ^— + V^\s) (126) 

s - Sa 

being k a constant. The above equation with the obvious replacements is identical to 
Eq. ([651) . As a final remark to this method, we add that the solution to the dispersion 
relation is not unique (any pair of solutions with different values of the constant k satisfy 
the same dispersion relation), since we may include any zeros to the inverse amplitude, 
which become poles in the amplitude. In the particular case of vrvr scattering such poles are 
inadmissible, so an additional condition on the constant k has to be imposed to describe the 
physics. 



C. Blackenbekler Sugar type equations 



It is interesting to see that many of our considerations are strongly related to a Blackenbekler- 
Sugar type equation [KH] . Let us consider the BSE in operator form. 



T = V + VGoT 



(127) 



As we know the "potential" V is the two particle irreducible amputated Green function. 
Thus, it does not contribute to the s-channel unitarity cut, i.e. Disc[V^fc(s)] = for s > 
4m^. On the other hand, the s-channel two particle reducible diagrams do contribute to 
the unitarity cut, but the part of them which does not contribute to the discontinuity is 
not uniquely defined. The idea is to split the two particle propagator into two parts, one 
containing the elastic unitarity cut and the rest. This separation is ambiguous, since we 
can sum to the former any function with no discontinuity. A practical way to do this is as 
follows. We write 



2 2 



G,{s) = A{q^)A{q_) = G,{s) + g,{s) 

where we impose that 

Disc[Gg(s)] = Disc[G'g(s)] = {-2m)6 
A solution to this discontinuity equation i^ 

Gg{s) = -m- — — — — 

qz _ _|_ _g^4 _|_ 



1 1 
q_ — m 



s > Am 



(128) 



(129) 



(130) 



Notice the identity {-2TTi)'^5+[q\ - m^]5+[q^_ - m^] = -2tt^5{P ■ q)5{q^ - m? + j) 
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This, of course, is not the only choice; we could equally well add a polynomial to Gq{s) 
without changing the discontinuity. Proceeding in perturbation theory we can see that 

T = V + V{Go + g)T = t + tGoT (131) 

where t = V + Vgt = V + VgV + VgVgV + . . .. 
Thus we are finally led to a equation of the form 

TAP, k) = tAp, k)+rrj ^Tp(g, k) Jjf^'^'^^ ^ ^^ Mp, g) (132) 

This equation looks like the Blackenbecler- Sugar equation |Q but with the important differ- 
ence that instead of the potential Vp{p, k) we iterate the reduced amplitude tp{p, k) which at 
lowest order in 1//^ coincides with the original potential. This equation satisfies de unitarity 
condition of Eq. (|). 

If we set the on-shell conditions P ■ p = P ■ k = we have a closed three dimensional 
equation for the amplitudes 0. 



D. Lippmann-Schwinger Equation 

The Lippmann-Schwinger equation has been employed recently Ref. ||^ in successfully 
describing s— wave meson-meson scattering using a coupled channel formalism]^. This ap- 
proach is also related to ours, although the way they obtain some results is not fully un- 
problematic. The authors of that reference compute the g° integral first, asumming that the 
contour can be closed at infinity, which for regular potentials and amplitudes may be accept- 
able. Then they introduce a three-momentum cut-off, hence breaking Lorentz and chiral 
invariance, and obtaining a three dimensional two body equation in the CM frame. Though 
this seems convenient for practical calculations, the justification for choosing this particular 
frame to include the cut-off is doubtful. Similarly to us, they also find that off-shellness 
may be ignored but it is not clear why chiral symmetry seems to play no role, particularly 
because their cut-off breaks it. In addition, to justify why the off-shellness can be ignored, 
they stick just to the one loop approximation. In our case we fully rely on the divergence 
structure of the non-lineal sigma model, to all orders in the loop expansion, and on its chiral 
symmetric structure and make use of the chirally symmetric dimensional regularization to 
establish that the off-shellness can effectively be incorporated as a renormalization of the 
two particle irreducible amplitude. 

Finally, in their finite cut-off framework they get the renormalization constant Cqo and all 
higher order low-energy parameter contribution to this channel, as a specific function of the 
cut-off. That function should be the same for all isospin-angular momentum channels. This 
introduces undesired correlations between the higher order corrections in different channels 



We have the formula d'^q = \\J {P ■ qY — dq^ d{P-q) cPq, where q'^ = —I = — and P-q = 0. 
Actually, their equation is the BSE, but we comply with their own notation. 
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which violates the spirit of ChPT. As a matter of fact, to reproduce the p— resonance within 
their framework requires a cut-off about three orders of magnitude bigger than that used in 
the (T— channel. We have already illustrated this point in Sect. |V] after Eq. ( |116| ). There, 
we compute the C parameter for the special case of dimensional regularization, but the 
discussion is identical if a three dimensional cut-off regularization is considered (see Eqs. 
(A3)-(A8) in the second item of Ref. [|). 



E. Inverse Amplitude Method 

Eq. (^) provides an exact solution for the tttt elastic scattering amplitude, thus it is 
not fully surprising that the lAM of Ref. 0] can be rederived from it. The essential point 
is our lack of knowledge on the two particle irreducible amplitude, V/-j, and the type of 
expansion proposed for it. For simplicity, let us consider the p— channel (/ = J = 1). In this 
channel, we showed in Subsect. |1V E| that if the Cu parameter is set to zero, the potential Vn 
has to diverge before the resonance energy is reached. That makes hopeless a perturbative 
expansion for V^i, but it suggests an expansion for 1/Vii, which should be small. Thus, 
expanding l/V/j (we come back to the general case IJ) in powers of 1/ P in Eq. ([7BD, we 
get 

Tu\s)\iAM = nT\'){s)]-' - (.)[t(;)(.)]-2 + ■ ■ ■ = -Us) + (133) 

where the Cu constant cancels out to all orders and V{sf/j^ = Tf] /{p — Tj^/tI'^). If t)j 
has a single non-kinematical Adler zero, like for s-wave scattering, then this zero becomes a 
double one at the first order of the approximation, since T = (T*^^))^/ i^pT^"^^ —T^^^). In other 
words, has a double pole, in contradiction with the dispersion relation for the inverse 
scattering amplitude, Eq. ( |123| ), as pointed out in Ref. 0. Notice that in the BSE the order 
of the Adler zero is always preserved. If in a particular channel, a resonance does not exist, 
one should expect 1/V not to be necessarily small and therefore the lAM might present 
some limitations. In our BSE framework, a chiral expansion of V is performed. Thanks to 
the inclusion of the renormalization parameter C, V might remain reasonably small in a 
wide region of energies. Besides, this scenario also allows us to compute form-factors. 



F. Fade Approximants 

In the Fade method 0, the series expansion in 1/p of the amplitude, 

TiAs) = W{s)/f + t\']{s)/P + Tf]{s)/f + ... (134) 
truncated at some finite order, 2N , is rewritten as a rational representation. 



r„W«rJf--'-'w = ^^i^ (135) 



where = 0, . . . , A^, Pf^-^^{s) and Q|f (s) are polynomials in 1/p of degrees 2N - 2K 
and 2K respectively. This approximation is enforced to reproduce the series expansion 
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for the amplitude and also to satisfy the unitarity condition, Eq. (pq). For the particular 
case 2N = 4, the only acceptable Fade approximant is the [2,2], and the lAM method is 
recovered. As we have already mentioned single non-kinematical Adler zeros are transformed 
at this level of approximation into double ones. In the language of the BSE, the Fade method 
is translated as a Fade approximant for the potential, for instance the [2, 2] approximant 
would be 

vir = ^f^^^^fL (136) 

which thanks to Eqs. (|65| ) and ( [74|) exactly yields to Eq. (|133|) . 

G. N/D Method 

In the N/D method one starts with the dispersion relation for the partial wave scattering 
amplitude (dropping the labels J, J for simplicity), 

r(3) = BM + ird/^''>'''>^-'''> (137) 

TT J4m2 S - S' 



where p{s) was defined in Eq. ( |12(]| ) and Disc[i?(s)] = for s > 4m^. B{s) presents a discon- 
tinuity for s < 0. The method is based on the following assumptions: T{s) = N{s)D{s)^^ 
and the discontinuity conditions 

Disc[Ar(s)] = D{s) Disc[T(s)] = D{s) Disc[5(s)] s < 

Disc[D(s)] = N{s) Disc[T"^(s)] = -N{s)2ip{s) s > Am^ (138) 

Dispersion relations for -D(s) and A^(s) with suitable subtractions determine the full ampli- 
tude solely in terms of the left hand cut discontinuity. The point here is that the disconti- 
nuity conditions in Eq. (|138|) do not determine uniquely the amplitude, as it also happened 
in the approach presented in Subsect. |VlIj| . Thus, one should supplement the discontinu- 
ity conditions with further information if possible. For instance the existence and position 
of the CDD (Castillejo, Dalitz and Dyson) poles as recently invoked in Ref. for 



meson-meson scattering, or of the bound states might help to find a physical solution from 
the whole family of solutions which satisfy the above discontinuity conditions. In any case, 
however, the multiplicative structure of the left and right hand cuts implied by the N/D 
method seems to contradict the additive cut structure deduced at two loops A 
more detailed study of the advantages and limitations of a N/D approach in this context 



will be presented elsewhere [57|. We note here, that the BSE preserves that additive left 



and right hand cuts, and hence a direct comparison with standard ChFT becomes possible. 

VII. CONCLUSIONS 

In this paper we have studied the consequences of chiral symmetry for tttt scattering and 
also for the scalar and vector form factors within an approach based on the BSE. Besides the 
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automatic incorporation of elastic unitarity, it becomes clear what is the subset of diagrams 
which is summed up. This requires the identification of a potential which corresponds to the 
amputated two particle irreducible Green function. As such, the identification is not unique, 
leaving room for a scale ambiguity, or equivalently a subtraction constant in the language of 
dispersion relations. We have dealt with the issue of renormalizability of the BSE within the 
framework of ChPT. We have recognized two alternative extreme viewpoints to implement 
the renormalization program which we have named: off-shell and on-shell schemes. 

In the off-shell scheme, the off-shell dependence of the potential is kept while solving the 
BSE. This way of proceeding allows to choose a renormalization scheme where the finite parts 
of the ultraviolet power divergences are not set to zero. This, still, produces more constants 
than those allowed by crossing at a given order in the chiral expansion, and it becomes 
sensible to fine tune them in a way to comply as best as we can with the crossing symmetry 
requirement. In practice we have found it convenient to implement this as a restriction on 
our partial wave amplitudes by matching those to Taylor expansions of the ChPT ones, 
including the contribution of the left hand cut. This off-shell scheme is rather simple at the 
lowest approximation level, and at the same time allows for a satisfactory fit to the phase 
shift TTTT scattering data in the / = 0, 1, 2 and J = 0, 1 channels leading to a prediction of 
the low energy parameters Ti, 2,3,4 with varying degree of accuracy [0. Conversely, using 
commonly accepted values for the /i,2,3,4 with their error-bars, a satisfactory and compatible 
prediction of the scattering data is achieved. In particular, our extrapolation of the low 
energy standard ChPT phase shifts in the I = J = 1 channel to higher energies predicts the 
correct p meson mass and width with 10% and 25% accuracy respectively. This approach, 
although very efficient to describe the data in the different isospin channels, becomes very 
cumbersome to pursue at next to leading order. 

Thus, we have also considered the on-shell scheme. As discussed in the paper, the fact 
that we are dealing with an EFT provides interesting insight into the problem. Given the 
infinite number of counter-terms necessary to ensure the renormalizability of the theory in 
a broad sense, one can reorganize the summation of diagrams in a way as to include in the 
potential the unknown coefficients. As an important consequence, a dramatic simplification 
in the solution of the BSE arises: all the divergences originated by the off-shell behavior of 
the potential can be effectively renormalized by redefining the potential and simultaneously 
ignoring the off-shellness in the solution of the BSE. In particular, on shell information 
for the full scattering amplitude can simply be obtained from the on-shell potential with 
unknown coefficients. We have found the on-shell scheme more convenient to take advantage 
of the ChPT information to one and two loops than the off-shell scheme. Eventually, both 
schemes would become equivalent if the exact amplitude were considered. 

In particular, by using the on-shell scheme we have constructed a tttt amplitude, which 
exactly reproduces the tree plus one loop scattering amplitude of ChPT and embodies exact 
elastic unitarity to all orders in the chiral expansion parameter 1//^. Thus, the definition of 
the phase shifts is unambiguous. In addition, a prediction for some ChPT two-loop low en- 
ergy parameters can be also made. The BSE within the on-shell scheme proves particularly 
useful when dealing with the vector and scalar form factors. Watson's theorem is auto- 
matically satisfied without explicitly invoking an Omnes representation. The requirement 
of an asymptotic behavior compatible with a once subtracted dispersion relation yields to 
a unambiguous prediction for the form factors, with no more undetermined constants than 
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those suggested by standard ChPT. Equipped with all this formalism we have been able 
to make a very accurate determination, in terms of estimated errors, of some low energy 
parameters. In particular, a remarkably accurate prediction has been achieved for the dif- 
ference li — I2 and the parameter Iq by fitting the vector form factor, rather than the vrvr 
scattering amplitude in the vector-isovector channel. Predictions for the corresponding 5ii 
phase shift with the propagated errors are satisfactory and accurate from threshold up to 
1200 MeV. Similar features are also found for other channels, although there the maximum 
energy, for which the predictions with their error-bars apply, ranges from 600 MeV in the 
scalar-isoscalar, to 1400 MeV in the scalar-isotensor channel, d— waves corresponding to the 
/ = and 2 channels are also describable up to 1000 MeV. In all cases we observe that 
the higher the energy the larger the error bars. This is not surprising since, our unitary 
amplitudes are generated from low energy information. 

From the point of view of predictive power we have also undertaken a thorough analysis of 
the predictions for the (9(1//^) contributions of the amplitudes to the threshold parameters. 



as compared to standard ChPT two-loop ones under a resonance saturation hypothesis [35 



The level of prediction is never less accurate than in the case of ChPT, and in some cases 
it is much better. The total predictions for the effective range parameters are in agreement 
with the experimental values within errors. 

The present calculation can be improved and extended in several ways which we describe 
in the following. First of all, the formalism presented in this paper can be enlarged to 
include coupled channel contributions. All of our formulae are valid for a coupled channel 
scenario just bearing in mind that now one is dealing with matrices instead of commuting 
C— numbers. This topic has been partially addressed in Sect. [IV C] . Secondly, our results, 
regarding the accurate description of one loop parameters with the present unitarization 
method, advises to implement the next order in our expansion. In this way we might try to 
determine, from fits to the data, more accurately some of the two loop parameters &i,2, 3,4,5,6) 
disposing the physically compelling, but unmotivated from a standard ChPT viewpoint, 
resonance saturation hypothesis. In this regard, the consideration of our formalism for 
describing Ki^ decays, in conjunction with the electric pion form factor might prove very 
fruitful in order to provide more accurate constraints for the one loop parameters. In neither 
case, have we addressed the determination of /s, which is related to the pion polarizability, 
and would require a thorough analysis of Compton scattering on the pion. In this respect, 
there exist already attempts in the literature trying to describe the physical process, but 
never focused from the point of view of determining this parameter. We leave these points 
for future research. 
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APPENDIX A: DETAILS ON THE OFF-SHELL BSE SCHEME 



1. 1=0 vrvr scattering 

The ansatz of Eq. ([TTD reduces the BSE to the hnear algebraic system of equations 

2/4(5) 



A 
B 
B 
C 



2/2 2/2 ' 2/2 

1 Ah{s)+Bh{s) 
P P 

1 (5m2 - 3s)/o(s) - 2/2(5) ^ (5m2 - 3s)/2(s) - 2/4(5 



B 



BI,{s) + Ch{s) 
P 



2/2 



B + 



2/2 



-C 



(Al) 



as we see there are four equations and three unknowns. For the system to be compatible 
necessarily one equation ought to be linearly dependent. This point can be verified by direct 
solution of the system. The integrals appearing in the previous system are of the form 



hnis) 



2\n 



(27r)4 [g2 - m2 + ig] [g2 _ ^2 ^ j^] 



(A2) 



/o(s),/2(s) and his) are logarithmically, quadratically and quartically ultraviolet divergent 
integrals. Translational and Lorentz invariance relate the integrals his) and his) with /o(s) 
and the divergent constants hi^^m"^) and /4(4m2) . 



his) = /o(4m2) + his) 



his) = ( 
his) = ( 



s/A)his) + hi^m^) 
s/AYhis) + hi4:m^ 



(A3) 



Note also that his) is only logarithmically divergent and it only requires one subtraction, 
i.e., his) = his) — /o(4m^) is finite ant it is given by 



his) 



1 



(47r) 



/ 4m2 
1 log 



4m^ 2 



(A4) 



where the complex phase of the argument of the log is taken in the interval [— vr, 7r[. Using 
the relations of Eq. ( [A3|) in the solution of the linear system of Eq. ([Al|) we get 



Ais) 


1 


■5m 


- Dis) 




Bis) 


-1 
-Dis) 


r 1 

[p 


Cis) 


1 


his 


Dis) 


P 



+ 



p p 

/2(4m2) + (m2 - s/4)/o(s) 



+ 



P 



(A5) 
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where, 



D{s) 



1 + 



/2(4m^ 



+ Io{s){ 



2s -m^ {s- Am^)l2{Am^) + 2h{Am^) 



2/2 



2f 



(A6) 



Eqs. (|A5|) and (|A6|) require renormalization, we will address this issue in Subsect. |III D 



2. 1=1 TTTT scattering 

The full off-shell scattering amplitude solution of the BSE in this channel is given by the 
ansatz of Eq. (^) with the functions M and given by 



2 / 2/;(4m') + (4m' - s)h(s) 
= 



-1 



sN{s) 



M{s) 4/2(4m2) - (4m2 - s)/o(s) 



6/2 l-/2(4m2)//2 
To obtain Eq. ( |A7| ) we have used that 



(2vr) 

1 pfJ. pv 1 

(^M^ ) [K - ^/4)/o(^) - H^ni^)] + -^?'^^/2(4m2) 

3 ^ s L -12 



(AT) 



(A8) 



3. 1=2 TTTT scattering 

The functions and C entering in the amplitude of Eq. ( p3D can be determined by 
solving the BSE in this channel, and thus we find 



where 



D{s) 



A{s) 
B{s) 
C{s) 

/2(4m2) 



1 rm2 /4(4m2) + (m2 - s/4)2Jo(s) 
1 ^ /2(4m2) + (m2 - s/4)Jo(s) 



1 lojs 
Dis) r 



P 



2m2-s (s-2m2)/2(4m2) +2/4(4m2 



2/2 



(A9) 



(AlO) 



Once again Eqs. ( |A9| ) and ( |A1(]| ) require renormalization, we will get back to this point in 
Subsect. pTD| . 
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4. On-shell and ofF-shell unitarity 



OfF-shell unitarity (Eq. (|^)) imposes a series of conditions to be satisfied by the amplitudes 
obtained in the previous subsections. Those read for the 1 = 1 case 



is — Arn^) 



12 

Disc [N{s)] 



Disc [N{s)] 
Disc [M(s)] 
and for the / = 0, 2 cases, they are 

Disc [A{s)] = - A{s) + B{s){m'^ - s/4) ^Disc [Io{s)] 



|iV(s)|2Disc [Jo(s) 



(All) 



Disc [C{s)] 
Disc [B{s)] 



B{s) + C{s){m^ - s/A) Disc[/o(s)] 



W - s/A)B{s) + A{s] 



m 



s/4)C(s) + B{s)) Disc [Io{s)] (A12) 



where Disc [f{s)] = f{s + ie) — f{s — ie), s > 4m^. 

On the other hand and thanks to the Schwartz's Reflex ion Principle, Disc [/o('S)] is given 

by 

Disc [/qI-s)] = h{s + i^) ~ h{s — ie) = 2iIm/o(s + ie) 



-i(27r) 



(2vr) 



'l _ l!!^, 5 > (A13) 



5. Renormalization method, crossing symmetry and lagrangian counter-terms 

The method of subtraction integrals makes any amplitude finite, by definition. In custom- 
ary renormalizable theories or EFT's, one can prove in perturbation theory (where crossing 
is preserved order by order) that this method has a counter-term interpretation. This is 
traditionally considered a test for a local theory, from which microcausality follows. Since 
we are violating crossing it is not clear whether or not our renormalization method admits 
a Lagrangian interpretation beyond the actual level of approximation. These are in fact 
a sort of integrability conditions; the renormalized amplitudes should be indeed functional 
derivatives of the renormalized Lagrangian. It is actually very simple to see that these con- 
ditions are violated by our solution. In terms of the generating functional Z[J] the four 
point renormalized Green's function is defined as 



1 S^ZlJl 

{0\T{7Ta{Xi)'!Tb{x2)'ITc{x3)'n:d{x4)}\0) - 



Z[J] 6Ja{Xi)6Jb{x2)6Jc{x3)6Jd{x4) 
MJc A^n MP r\'^P' 



J=0 



(27r)4 (27r)4 {2%)^ {2tt)^ 
X Atb'ip-)i-i)Tp{p, k)a'b';c'd'5\P - P')/\,,.{k+)AM>{k^) (A14) 



54 



The integrability conditions are simply the equivalence between crossed derivatives which, 
as one clearly sees, implies in particular the crossing condition for the scattering ampli- 
tude. Since crossing is violated, the integrability conditions are not fulfilled and hence our 
renormalized amplitude does not derive from a renormalized Lagrangian. 

Crossing symmetry can be restored, in certain approximation, as it is discussed in Sub- 
sect. pTE| , by imposing suitable constraints between the subtraction constants which appear 
in the renormalization scheme described in Subsect. |111 D| . These conditions reads: 

75J2^'^=V2m^ + 8/0^'^=^ + 33/0^'^=° + 5/,^'^=Vm^ + = 

rRj=2 8Jo^''=° 719 _^ 



15 5 72007r 



2 



4m2 25 12m2 360007r2 

171^''=' 111/0^'^=° 3/,^'^=^ _ 2159 ^ 

16m4 500 4000 200m2 40m4 4800007r2 ^ ^ 



The first of these relations was obtained in Ref. 12 



Finally and for the sake of clarity we quote here the result of inverting the Eq. (18) of 
Ref. [g. 



/o''''"' = -Y^ (2(^2-^1) +97/60 



2 

m 



^2''-' = ^ ((2(^2 - k) + - 65/24) 

lo''^' = {22h + 2812 + 31/2) 

4 

jR,i=o ^ m f^^2h - 56812 + 60013 - 672r4 + 1057) (A16) 
76807r^ ^ ^ 



6. Remarks on the choice of the pion field 

The solution of the BSE requires an ansatz for the off-shell potential Vp{p, k). The non- 
linear realization of chiral symmetry through the non-linear a-model posses an additional 
difficulty. The pion field, 0, is encoded in the unitary space-time matrix U{x), U^U = 1, 
so it lives in the three sphere of radius /. However, the particular set of coordinates is 
not unique. For instance, one may have a polar representation U{x) = e"^"^/-^, being r the 
Pauli matrices , a cariesmn representation U{x) = — (p'^/p + icf-r/f, or a stereographic 
representation U{x) = (l + i0-r/2/)/(l— i0-r/2/). It is well known that any representation 
yields to the same kinetic and mass terms, but different pion interaction terms. Indeed, up 
to second order in the chiral expansion, the effective lagrangian reads 

i'2 £2 2 

C = LtT{d^U^d^u) + ^-^tr(u + f/t - 2) (A17) 
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which up to fourth order in the pion field becomes 



C = lidjr -^P- ^,idJrP+ ^-^{^djf + {a- (A18) 

where the unfixed a parameter has its origin in the arbitrariness on the form of the ?7— matrix. 
Thus, the polar and cartesian representations correspond to a = whereas the stereographic 
representation leads to a = 1. The a dependence disappears if one uses the equation of mo- 
tion. The above Lagrangian gives 

Ap{p, k) = --1 |s (^1 - I) - a (f + e) + {2a - 1)| (A19) 

which by means of the symmetry properties discussed in Sect. |1| completely defines the 
two identical isovector (off- or on-shell) meson scattering amplitude {Tp{p, k)ab;cd) for the 
process {P/2 + p, a) + {P/2-p, b) {P/2 + k,c) + {P/2 -k,d), defined in Eq. (|). On the 
mass-shell the dependence on a of the amplitude disappears. Off the mass shell an explicit a 
dependence is exhibited. Thus, the BSE kernel Vp{p, k) will depend on a and might generate 
also an a— dependence in the solutions of the BSE. For on-shell pion scattering and when 
the exact potential V and propagator A are used the a dependence drops out. However, 
when approximated V and A are inserted in the BSE, the solution of the equation might 
display an a dependence, even for on-shell pions. This is the case for the lowest order / = 
and 1 = 2 on-shell solutions, 

T-'(s) = -I (s) + 2{f + {l-5a/2)h{Am')f 

°° (m2-2s)/2 + (l-5a/2)2(2/4(4m2) + (s-4m2)/2(4m2)) ^ ' 

T-^(,) = r , 2{p + {l-a)I,{4m^) f 

2° ^ ^ ^ (s - 2m2)/2 + (1 - a)2 (2/4(4m2) + {s - Am^)l2{4m^)) ^ ' 



where both amplitudes above require renormahzation and reduce to those given in Eqs. 
and Eqs. (|2^) for a = 0. The lowest order J = 1 BSE solution turns out to be independent 
of a. 

It is worth noticing that the energy dependence encoded in the amplitudes above is not 
changed by this a dependence, and they can generically be written as 

7^0 (20) (^) = -^o(^) - c + (A22) 

with a, h and c suitable constants. So there is redundancy of parameters, since we have 
four parameters Jo(4m2), J2(4m2), J4(4m2) and a, mapped into three parameters a, 6 and 
c. Obviously, the particular value of a would be irrelevant, when a fit of these constants to 



data is performed, as it was done in Ref. [12|. The trouble might appear, as we will see. 



when matching these amplitudes to those deduced in ChPT, which of course do not present 
any dependence on a for pions on the mass shell. 

The a dependence shown in Eqs. ( |A20|) and ( A21 ) is undesirable because we are dealing 



to on-shell amplitudes. It appears because we have not generated loops in the t— and 
M— channels, i.e. the amplitudes do not satisfy exact crossing. On the light of this discussion, 
we envisage two alternative procedures linked to different renormahzation schemes: 
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Within the renormahzation scheme presented in Sect. |IIID| , we can match the 1//^ 
expansion of the BSE amphtude up to (9(1//^), to the threshold expansion, up to 



order (s — 4m^)^, of the Gasser-Leutwyler amphtudes, as it has been done in Sect. |A 5 



This procedure modifies the last of Eqs. (|A16|) and all of the constraints of Eq. ([A 151) 



and ensures the a independence of the matched terms, but there remains a residual a 
dependence in the renormalized on-shell BSE amplitude, starting at order 0{l/f^). 
Actually, this matching procedure does not work ifa = 2ora = — 1, because in this 
cases the BSE coefficients of {s—Aim?)/f^ and {s—Am^Y/f^ pieces are not independent 
for / = and 1 = 2 respectively, making then the matching to the Gasser-Leutwyler 
amplitudes overdeterminecQ. In addition, we do not have any restriction on the 
possible a values. This is way we find more appropriate the following scheme. 



From the renormahzation discussion presented in Sect. [111D| it is clear, as it is also 



the case in standard ChPT, that to achieve renormahzation of the amplitudes an 
infinite number of counter-terms is required. Being then the finite parts, order by 
order in the perturbative expansion, undetermined and have to be either fitted to data 
or determined from the underlying QCD dynamics. This implies in particular that 
the divergent integrals appearing in the numerator and denominator in Eqs. ( |A20|) 



and ([A21|) may be chosen to be independent of each other. This can be used ad- 



vantageously to eliminate the explicit a dependence, by simply redefining the inte- 
grals, /2(4m^)(l — 5a/2) — > /2(4m^), /2(4m^)(l — a) ^ /2(4m^) in the numerators 
of Eqs. ( [X20D and ( [MID , and /2,4(4m2)(l - 5a/2)2 hA^m^), /2(4m2)(l - a)^ 
/2(4m^) in the denominators of Eqs. (|A20| ) and ( |A21|) , respectively, since we know that 
the total amplitude must be a independent. Even though this choice might appear 
more arbitrary than the scheme presented in Sect. PIIE)| , it is still more restrictive than 
what the renormahzation of an EFT allows for. 

From the discussion above, we prefer to work within the second scenario where the 
unphysical a— dependence can be ignored, when looking at on-shell amplitudes, and the 
particular value a = can be considered, as we have done along Sect. |T|. Within this 
renormahzation scheme any value of a would lead to the same results presented in that 
section. 



APPENDIX B: LEADING AND NEXT-TO-LEADING ELASTIC tttt 
SCATTERING AMPLITUDES IN CHPT 

The O(p^) + 0{p^) on-shell SU{2)— ChPT elastic tttt amplitudes, expressed in terms of 
the renormahzation invariant parameters /i, 2,3,4, are given^ by 



29That would require a restriction among the Vs. 



^°Note that the function /o(s), fulfilling loiim^) = is related to that of (T|] fulfilling J(0) = 
by 
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Tij{s)^Trj{s)lf + Trj\s)/f 
t\'}{s)^t\']{s) + Us)^[T?}{s)] 

I ^0; J^O 



(Bl) 
(B2) 



m'^-2s 







6 







1 = 
7 = 1 
1 = 2 
1 = 2 



J = 2 
J=l 
J = 
J = 2 



(B3) 



1 1 , 



(B4) 



hij = - 



1 



2 7- 



5?Ti^ I 101m^(s-4m^) , 191(s-4m^)^ , 
47r2 967r2 2887r2 

287(5-4m^)^ I 481(5-4m^)^ 
28807r2 67200m2 7r2 



d(cos^) {fi(t)Io(t) + fi(u)Io{u)) Pj(cose) 

7 = 0; J = 

7 = 0; J = 2 

7= 1; J= 1 

llm^ I 179m^(s-4m^) , 769(s-4m^)^ , r _ o • 7" — H 

47r2 967r2 14407r2 1 — Z , J — U 

529(a-4m^)^ , 277(5-4m2)3 , T _ n . j _ r, 

144007r2 ^ 67200m2 7r2 J — Z , J — Z 

where i = -2(f - m2)(l - cos^), w = -2(f - m2)(l + cos^) and 



89m2(s-4m2) _ 37{s- 
— < 2887r2 288 



28807r2 + ■ ■ ■ 



(B5) 



' lOx^ + x{2s - 32m2) + 37m^ - 8sm^ 7 = 
2x'^ + x{s + 2m^) — m"^ — 4sm^ 7 = 1 

^ - x(s + 2m2) - 5m^ + Asm'^ 1 = 2 



(B6) 



7o(s) = -J(s) + ^ 
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' (40/1 + 8O/2 - 15/3 + 84/4 + 125) + m2(s - 4m2) x 



gij 



4m2 (-2/1 + 2/2 + 3/4 + 1) + (s - 4m2) (-2h + 2/2 + 1 



{l6h + 32/2 - 6/3 - 24/4 + 50) + m^{s - 4m2) x 



/ = 0; J = 
/ = 0; J = 2 
/ = 1; J= 1 

/ = 2; J = 
I = 2; J = 2 



(B7) 



The integrals which appear in the definition of hjj in Eq. ( |B4|) can be done analytically 
by means of the change of variables: 



cos 6 —>■ X 



1 4m'^ I -1 



- 1 



1 



C(cose) 



(B8) 



being ^(cos6') = t or u. Thus we get, 

/i/j(s) = ajj{s) + bij{s)L{s) + cij{s)L{sy 

where 



L{s) = log( 



1 + ^/l _ 



1 - a/1 - 



4m^ 



(B9) 



(BIO) 



and s > Am?. In the following we list some special cases of interest, 



&00 
Coo 

an 
fell 

Cll 



-506m^ + 130m2s - lls^ 
1447r2 

75m^ - 40m2s + I 



s — Am? 



247r2 
m^(25m2 — 6s) 
87r2(s-4m4) 

120m^ - 149mS + 37m2s2 + 

1447r2(s - 4m2) 
36m*^ - 72m^s + IGm^s^ - T 
487r2(s - 4m2) \] 

m^{6m^ + 13m^s — 3s^) 
87r2(s - 4m2)2 



(Bll) 



4m2 



(B12) 
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308m^ - 58m^s - 25s'^ 

2887r2 ' 

6m^ - 32m2s + lls^ I s 

m^iw? + 3s) 
87r2(s - 4m2) 

-31936mS + 27646m6s - 4371m^s2 + Qlw?s^ + 14s^ 

7207r2(s - 4m2)2 
516m® - 760m'3s + 180m^s2 _ 13^2^3 ^ ^.1 i - 

247r2(s - 4m2)2 s - Am? 

m^(172m^ - 98m^s + 49m2s2 _ 65^) 
87r2(s - 4m2)3 

124672m® - 48832mS + 53592mS2 - 5182m2s3 + 247s^ 

144007r2(s - 4m2)2 
-480m® - 980m6s + 117mS2 - 36m2s3 + 2s^ / i 

1207r2(s - 4m2)2 \ s - Am? 

m^(32m6 - 76mS + Ilm2s2 _ 3^3) 
87r2(s - 4m2)3 
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